NASA TN D-8151 


NASA TECHNICAL NOTE 



NASA TN D-8151 


FUNDAMENTALS OF FLUID SEALING 


John Zuk 

Lewis Research CeMer 
Cleveland f Ohio 441 3 5 





NATIONAL AERONAUTICS A ND SPA CE ADMINISTRATION • WASHINGTON, D. C. • MARCH 1976 






1 . Report No. 

NASA TND-8151 


2. Government Accession No. 


4. Title and Subtitle 

FUNDAMENTALS OF FLUID SEALING 


3. Recipient's Catalog No. 


5. Report Date 

March 1976 


6. Performing Organization Code 


7. Author(s) 

John Zuk 


8. Performing Organization Report No. 

E-6910 


9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 


10. Work Unit No. 

505-04 


11. Contract or Grant No. 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 


13. Type of Report and Period Covered 

Technical Note 


14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

The fundamentals of fluid sealing, including seal operating regimes, are discussed. The general 
flu id- flow equations for fluid sealing are developed. Seal performance parameters such as 
leakage and power loss are presented. Included in the discussion are the effects of geometry, 
surface deformations, rotation, and both laminar and turbulent flows. The concept of pressure 
balancing is presented, as are differences between liquid and gas sealing. Also discussed are 
mechanisms of seal surface separation, fundamental friction and wear concepts applicable to 
seals, seal materials, and pressure-velocity (PV) criteria. 


17. Key Words (Suggested by Author(sM 

Seal; Face seal; Gas-film seals^ Labyrinth 
seals; Lubrication; Incompressible flow; 
Compressible flow; Narrow slots; Noncontact- 
ing seals; Fluid- film seal^ Hydrodynamic 
seals; Hydrostatic seals; Friction; Wear; 
Boundary lubrication; Seal materials 

18. Distribution Statement 

Unclassified - unlimited 
STAR Category 37 (rev. ) 

19. Security Classif. (of this report) 
Unclassified 

20. Security Classif. (of this page) 
Unclassified 

21. No. of Pages 
166 

22. Price* 

$6. 25 


For sale by the National Technical Information Service, Springfield, Virginia 22161 












CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

SEAL OPERATING REGIMES 2 

GENERAL EQUATIONS OF MOTION FOR VISCOUS FLUID FLOW 4 

PHYSICAL QUANTITIES OF INTEREST 9 

Leakage Flow Rate 9 

Pressure Distribution 10 

Opening Force • • 19 

Frictional Horsepower Requirements 10 

Film Stiffness H 

Center of H 

Film Temperature Rise 11 

SEAL PRESSURE -BALANCING FUNDAMENTALS 12 

LAMINAR FLOW 14 

Flow Between Parallel Plates in Relative Motion 14 

Flow Between Flat Convergir^ or Diverging Plates 20 

Axial Flow Between Annular Cylinders 27 

TURBULENT FLOW 29 

Turbulent Flow Between Parallel Plates 31 

Turbulent Flow with Small Linear Deformation (and Constant Width) 32 

Turbulent Radial Flow 32 

Turbulent Flow Between Concentric and Eccentric Cylinders 33 

APPLICATION OF THEORY TO FACE -SEAL CONCEPTS 34 

COMPRESSIBLE FLOW OF GASES 34 

Compressible Flow Equations ...... . . . 35 

Classical Viscous Compressible Flow Model 36 

Pressure Distribution - Application Example 43 

Approximate Models of Compressible Flow 43 

Geometry Effect on Leakage 48 

4Q 

Openii^ Force 

Center of Pressure 49 

Continuvuu and Noncontinuum Flow 50 


iii 


Flow Throt^h Porous Media 51 

INVISCBD FLOW EQUATIONS 51 

Gravitational-Head Orifice Flow 53 

Nozzle Flow 53 

Flow Function Approach 55 

Labyrinth Seals 56 

ENTRANCE FLOWS AND LOSSES 57 

FLOW REGIMES 60 

MECHANISM OF FILM PRESSURE GENERATION BETWEEN RELATIVELY 

MOVING SURFACES 64 

Hydrodynamic Pressure Development 64 

Squeeze-Film Concept 68 

Cavitation and Degasification 74 

Pressure Development due to Microroughness 75 

Pressure Generation due to ’’Wavy” Surfaces 76 

BOUNDARY LUBRICATING REGIMES 76 

Fimdamental Friction Concepts 76 

Fundamental Wear Concepts 78 

SEAL MATERIALS 82 

CONCLUDING REMARKS 83 

APPENDIXES 

A - SYMBOLS 85 

B - VISCOSITY OF FLUIDS 90 

C - COMPRESSIBLE FLOW - AN ALTERNATE THEORETICAL APPROACH TO 

FINDING THE ORIFICE EXPANSION FUNCTION 92 

D - DERIVATION OF REYNOLDS LUBRICATION EQUATION 94 

REFERENCES . 98 


iv 



FUNDA/VIENTALS OF FLUID SEALING 
by John Zuk 
Lewis Research Center 


SUMMARY 

The fundamentals of fluid sealing, including seal operating regimes, are discussed. 
The general fluid-flow equations for fluid sealing are developed. Seal performance pa- 
rameters such as leakage and power loss are presented. Included in the discussion are 
the effects of geometry, surface deformations, rotation, and both laminar and turbulent 
flows. The concept of pressure balancing is presented, as are differences between 
liquid and gas sealing. Also discussed are mechanisms of seal surface separation, fun- 
damental friction and wear concepts applicable to seals, seal materials, and pressure- 
velocity (PV) criteria. 


INTRODUCTION 

In fluid sealing as in every branch of engineering, one must draw on the experience 
and knowledge of many other disciplines and fields. The ei^neerii^ of seals can involve 
fluid mechanics, heat transfer, lubrication theory, structural and solid mechanics, 
thermodynamics, chemistry, physics, metallurgy, and dynamics, as well as other 
fields. Seal problems may consist of a superposition of effects which can be interre- 
lated. Usually each effect can be analyzed by itself. Then integrated effects must be 
evaluated. 

Fluid sealing is generally the same as any other branch of engineering except for 
the importance of small-scale effects. Seals are characterized by surfaces in relative 
motion separated by a very narrow gap. In order to maintain proper operation, very 
small differences in the dimensions of seal parts must be maintained. Deformations in 
geometry due to imposed thermal gradients, frictional heating, pressure gradients, and 
mechanical and inertial forces must be held to a minimum. Actually, in some cases, 
deformations must be no more than microvalues. 

The fundamentals of fluid flow are important to understanding the various sealing 
devices. This presentation discusses those fundamentals as they apply to seals. Basic 



principles of incompressible (liquid) and compressible (gas) sealii^ flows are studied. 
The mechanism of film pressure generation between relatively moving surfaces is de- 
scribed. Fundamental friction and wear concepts, various seal lubrication operating 
regimes, and surface topography effects are also presented. 

Generally, sealing flow is a direct opposite of conventional fluid flow through pipes, 
ducts, and nozzles, where it is desired to have as efficient a flow process as possible 
(low friction). In seals it can be desirable to have an inefficient flow process because 
this results in low leakage; however, a trade-off must be made with power loss. 

This presentation is intended to provide the backgroxmd necessary for a fimdamental 
understanding of fluid sealing in general and of specific seal types in particular. Unfor- 
timately, this discussion of the theory leaves much of the subject uncovered. Deriva- 
tions are omitted unless they are essential to the understanding or illustration of a con- 
cept. A remark, perhaps, should be made concerning the accuracy of the applicable 
equations contained herein. When we are dealii^ with fluids and seal configurations for 
which test data are available, design accuracy of perhaps 25 percent may be realized. 
When we are dealing with different configurations or using fluids that deviate from ideal 
properties, present knowledge is rated as fair or poor and testing becomes necessary in 
the design. 

Fundamentals of fluid sealing can also be found in a textbook by Mayer (ref. 1). A 
survey report (ref. 2) should be useful. Each topic covered includes, where possible, 
a reference wherein more information can be found. The seal nomenclature used is 
that adapted by the ASLE Seals Technical Committee (ref. 3). 

In all references cited, the U. S. customary system of units was used. The Inter- 
national System of Units (SI) has been added to the text and figures for this report only. 


SEAL OPERATING REGIMES 

Seals operate in many lubrication regimes depending on the type of seal, the sealed 
fluid, the application, and so forth. A plot of friction coefficient against seal duty pa- 
rameter illustrates the various seal operatii^ (lubricating) regimes that can exist 
(fig. !)■ 

For illustrative purposes, consider a lift-off type of seal that is in rubbing contact 
at startup and shutdown. The way this seal may change from one lubricating regime to 
another within an application is shown in figure 1 for liquid lubrication. The figure 
shows the variation in friction coefficient a seal imdergoes from startup under a load 
(e. g. , spring force and pressure) in the boundary lubricating regime to the steady-state 
operatii^ speed in the full-film lubricating regime. (The mechanism for achieving full- 
film operation could be an external pressurizaticm source, or it could be self-generated 
by hydrodynamic lubrication. At startup the parts are in solid- to- solid contact and the 
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seal runner (rotating member) begins to turn under essentially dry conditions and starts 
to follow the path AB. If sufficient lubricant is available, the luTsricant is ordinarily 
drawn between the sliding surfaces at once (by capillary action or a forced pumping 
mechanism) and the seal immediately enters the thin-film (or mixed film) regime, fol- 
lowing the path ABC. When the speed reaches the value corresponding to point C, the 
seal enters the full-film regime, in which it remains until comir^ to operating speed at 
point D. It can be seen in figure 1 that, if the lubricant were not present, the seal would 
be forced to operate dry at a speed corresponding to point A. The resulting temperature 
rise could be extreme because of the high friction. 

Now some of the details are examined more closely. The curves shown in figure 1 
were originally for journal bearings, but the same principles apply for a seal. Friction 
coefficient is plotted against a seal duty parameter juN/P^, where ii is the fluid vis- 
cosity, N is the rotational speed, and is the net seal-face load. (All symbols are 
defined in appendix A. ) To the right of the second dashed vertical line is the region of 
full-film fluid lubrication, that is, thick-film lubrication, where the surface asperities 
are completely separated by an oil film of such thickness that no metal-to-metal contact 
can occur (fig. 2(a)). The friction coefficient for a hydrodynamic film can be calculated 
from 


Net normal closir^ force 

where tA is the traction force. Hydrodynamic lubrication theory applies, and the flow 
is laminar. At sufficiently large values of the seal duty parameter, turbulent flow can 
occur (transition occurs at point D and turbulent flow exists beyond point E of fig. 1). 

The friction here rises significantly and increases at a more rapid rate with speed than 
in the laminar flow regime. To the left of the dashed vertical lines are the regions of 
boimdary and thin-film lubrication. As noted in figure 2(c), the film thickness in bound- 
ary lubrication is so small that asperities make contact through the fluid film. The thin- 
film regime is one that combines hydrodynamic and boundary lubrication. This is also 
the regime where elastohydrodynamic effects may be important. Friction coefficients 
in the boimdary and thin-film regimes are empirically determined. 

In full-film fluid lubrication, since the asperities do not contact, only bulk lubricant 
physical properties are important. In boundary and thin-film lubrication, the properties 
of the metals as well as surface physics and chemistry are of primary importance since 
there is metal-to-metal contact by asperities. Lubricant chemical properties can in- 
fluence the type of damage that occurs. 

The lubricating regimes can similarly be associated with seal-face loading and 
speed. The three regimes are shown from this point of view in figure 2. That is, the 
full-film lubricating regime is characterized by the film thickness being several times 
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greater than the surface roughness. The thin-film regime has the film thickness of the 
order of the surface roughness. In the boundary regime, asperity contact characterizes 
the interface. 

The friction coefficient behavior is different for a gas (fig. 3). Since gases are poor 
boundary lubricants, the friction coefficient values are almost those obtained for a solid 
sliding on a solid at low seal duty parameters. Since a gas has a much lower viscosity 
than a liquid, friction forces can be one or two orders of magnitude less in the full-film 
lubricating regime. For a gas film seal to operate in this regime, incorporation of a 
lift geometry to the seal faces is usually required. Operating gaps are inherently smal- 
ler for gas film seals (due to the low viscosity); therefore, these seals are more sensi- 
tive to face distortions. This, coupled with the poor boundary lubricating properties of 
gases, means that stable self-induced hydrod 5 mamic operation is unlikely in gas film 
seals. The friction coefficient variation is qualitatively similar to that for liquid lubri- 
cation in the full-film regime until compressibility effects become significant. Gener- 
ally, compressibility effects become important before turbulence occurs. However, for 
large gaps and sufficiently high pressures or speeds, turbulent flow may occur. 


GENERAL EQUATIONS OF MOTION FOR VISCOUS FLUID FLOW 

A brief derivation of the fluid-flow equation governing many seal flow situations will 
now be given. More details can be found in any fluid mechanics textbook, such as refer- 
ence 4. However, we will follow the development given in reference 5, which derived 
the equations for seal applications. 

Consider the elemental fluid particle dx dy dz moving in a pressure field, as 
shown in figure 4. The element dx dy dz is small enough that we can use differential 
and integral calculus but large enough for a representative statistical average. That is, 
the density of the fluid varies smoothly. The general conditions for equilibrium will be 
formulated; no stipulation will be made as to the nature of the fluid other than that it is 
viscous and follows Newton’s law of viscous flow. (Viscosity and both Newtonian and 
non-Newtonian fluids are discussed briefly in appendix B. ) Newton’s law of viscous flow 
states that the shear stress resulting from the shear deformation is a linear function of 
the rate of deformation. That is, the fluid property is a fvmction of the thermodynamic 
state only. Mathematically in rectilinear Cartesian coordinates, this can be expressed 
as 
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( 2 ) 



This is an empirical observation. It is the only empiricism used in deriving the equa- 
tions of motion. 

According to the conservation- of-momentum principle, which is Newton's second 
law, dynamic equilibrium of forces on the element dx dy dz requires that 
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Figure 5 illustrates the elemental volume forces acting in the x-direction: 
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( 6 ) 
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Each of the pressure terms in equations (4) to (6) are normal stresses and can be shown 
to be a fimction of an average pressure P and velocity gradients. The result is 


p =-p-2 2^1 h 

3 V3x By Bzl Bx 

p . -p-2 + 2^3«r 

y 3 \3x 9y 0z/ dy 

p =.p.2^/iu^3w^9v\ 2^ii 

^ 3 \0x 3y 9z/ az 


(7) 

( 8 ) 
(9) 


Equations (7) to (9) are the constitutive equations for a Newtonian fluid. The pressure P 
is the thermostatic (isotropic) pressure. If equations (1) to (3) and (7) to (9) are suc- 
cessively differential and substituted into equations (4) to (6), the following equations re- 
sult: 


where 
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is the material or Euler ian derivative. It is composed of the local acceleration (tem- 
poral velocity change) and the convective acceleration (spatial velocity chaise). The 
convective acceleration is due to the coordinate system, which is fixed relative to fluid 
motion (laboratory coordinate system). Thus, the Navier-Stokes equations are written 
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in the Eulerian coordinate system. 

Equations (10) to (12) are the Navier-Stokes equations, which are the most general 
equations of motion for viscous flow. For specific problems they are simplified, which 
hopefully results in equations that can be solved. In this general form they are much too 
complex to solve even on the most advanced computers. Solutions to specific problems 
require satisfaction of the boundary conditions and satisfaction of the continuity equation 
(conservation of mass). The continuity equation can similarly be derived as the momen- 
tum conservation equation and is 

dp ^ 3(pu) , a(pw) 3(pv) _ n (13) 

0t 3x 9y 9z 

This is the differential form of the mass balance equation for a control volume. 

The continuity equation with the three Navier-Stokes equations yields a well-posed 
mathematical problem, that is, four equations and four unknowns for liquids (appropri- 
ate boimdary conditions must be specified). 

Now the equations will be further simplified. The resulting form will be mathemat- 
ically tractable and yet will retain the important physics of sealii^ flows: 

(1) For steady flow. 


^ ^ _ 9v _ Q 

at at at at 

(2) Generally, the flows can be described as two dimensional; hence. 


3y 

(3) The flow is constrained to flow in one direction only, which can arbitrarily be 
assigned to be parallel to the x-direction; thus, v = 0 with its derivatives (fig. 6). The 
planes are considered to be infinite in the y-direction, but practically they are of finite 
width W. Any end-wall effects are thus neglected. The velocity must be zero at the 
boimdaries because of the physical "no slip” of the fluid at the boundary. Since the 
z-dimensions are much smaller than the x-dimensions, only the velocity derivatives 
taken with respect to z are appreciable. Therefore, 

^ ^ a^u\ ^ 0 
W’ ax2’ ay2; 
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In the general case of converging or diverging channels, v and its derivatives may be of 
appreciable magnitude. However, the assumption that they are vanishingly small has 
practical validity in seal performance analysis since converging or diverging angles of 
actual seal faces are, of necessity, extremely small. 

Assume that the fluid is incompressible, which is the case for liquids (not true for 
gases). That is. 


p = Constant 

The resulting simplification of equations (10) to (13) yields 
Momentum conservation; 


X--— +i^-— = 0 (14) 

P 9x p 3^2 

Z - - — = 0 (15) 

p 9z 

Y = 0 (16) 

Mass conservation: 

iH+iX = o (17) 

9x 9z 

We can assume that the body force Z is due only to the static (gravitational) head so 


Z = — (gz) = g 
9z 


( 18 ) 


Thus, 


1 9P _ 1 dP _ 

o 

p 9z p dz 


Hence, 
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P = / pg dz = pgh + f(x) 

*'o 


( 19 ) 



This indicates that the z-direction pressure variation is strictly hydrostatic, which is 
consistent with the unidirectional, two-dimensional flow assumption. 

Equations (14) to (17) are the classical forms of the incompressible, purely viscous 
flow equations. A number of revelant cases can be solved from this set. Li the section 
LAMINAR FLOW, details will be given for flow between parallel plates in relative mo- 
tion because this illustrates the basic concepts. Other important cases will follow with- 
out details. Details can be carried out readily as an exercise or by referring to the 
cited references. In the next section the usual physical quantities of interest will be 
described. 


PHYSICAL QUANTITIES OF INTEREST 

Typically, in fluid-film sealing the pressurized fluid conditions and ambient condi- 
tions are given (e. g. , sealed pressure and temperature). Geometric envelope con- 
straints are almost always known, although it is usually desirable to refine the geometry. 
The geometric orientation may also be dictated by system considerations (e. g. , radial 
face seal or circumferential seal). Information on seal leakage flow rate, pressure dis- 
tribution, separating (or opening) force, frictional horsepower requirements, film stiff- 
ness, center of pressure, and film temperature rise is usually desired and will be 
evaluated herein for most common cases. Figure 7 illustrates a pressure-balanced face 
seal which can be used as a reference. 


Leakage Flow Rate 

The desired information may be the volumetric flow rate Q for incompressible 
fluids, where 


or mass flow rate 


Q= /u 


dz dx 


M = pQ = pUg^yA 


This is the measure of seal leakage. 
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Pressure Distribution 


Pressure distribution is especially necessary for pressure-balanced seals. It is 
usually found from an analytical expression which is a function of the flow rate. 


Openii^ Force 

Since there is a pressure drop when a fluid flows, the integrated effect yields a net 
force. This force is called the separating force (also called opening force). This force 
can be found by integrating the pressure over the total surface that the pressure acts on. 
In seal work, it is generally desirable to know only the increment of lifting force due to 
the pressure differential rather than the absolute pressure force. Hence, the pressure 
force due to ambient (reference) conditions must be subtracted. Thus, 

fnel= I P'JA-P^A 


Frictional Horsepower Requirements 

Power requirements are also important. Seal friction can be generalized to the ex- 
tent that the surfaces in contact may be either ’’dry” or ”wet. ” In either case the fric- 
tional horsepower requirements are found from 

2ffNT 
~ 396 000 

In dry friction, where F is the net seal-face force acting between the two rubbii^ sur- 
faces at some effective hydraulic diameter D, the tangential force opposing rotation is 
the product of F and the coefficient of Coulomb friction tj. For wet friction, where a 
viscous film separates the surfaces, the horsepower requirement is predicated on the 
rate of energy dissipation in viscous shear. The common assumption of Newtonian fluid 
applies, such that 


3u 

T = jU — 
0Z 


Thus, the torque is 
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TX dA 


Film Stiffness 

The film stiffness is important for maintaining positive surface separation when the 
seal is operating under severe conditions, such as high sealed-fluid pressure differen- 
tials and/or speeds. The axial film stiffness is defined as 


dF 


K = 


net 


dz 


Center of Pressure 

The center of pressure is important in determining the net twisting moment that the 
fluid film exerts on the seal nosepiece. The material of the rotor must accommodate the 
residual stress due to this moment or deflect in a satisfactory way. The net twisting 
moment can be found from 


X 


c 


/(P - P^)x 


dx 


z. 

w 


Film Temperature Rise 

The film temperature rise due to viscous shearing is also important. It is desirable 
to know this value to assure that the seal can tolerate this temperature rise. A very 
rough estimate of the film temperature rise due to viscous shearing can be made by 
equating the heat generated by viscous shearing with the heat transferred by convection. 
Thus, 


T 

film, av 


- T 

seal surface 





CpM 


This calculated film temperature rise will be higher than the actual film temperature 
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rise. The predominant mode of heat transfer is usually conduction by the seal surfaces, 
not convection. 

Seal fluid flow can be analyzed by usii^ either an exact or an approximate analysis. 
Seal flows which are purely viscous laminar flows are amenable to a differential analy- 
sis. This type of analysis yields finely detailed behavior of the fluid flow. Unfortu- 
nately, there are situations where the differential analysis yields equations which are 
impossible or impractical to solve for design analysis purposes. Such cases include 
turbulent flow, where exact physical knowledge is unknown, and choked flow, where non- 
linear behavior characterizes the flow. In these cases, approximate methods, such as 
integrated average methods, must be used. Although the integral models satisfy only 
mean conditions in the flow field, they yield good results on gross quantities, such as 
seal leakage and pressure distribution. 


SEAL PRESSURE-BALANCING FUNDAMENTALS 

One of the primary objectives in fluid-film, face-seal design is to ensure that the 
face loadii^ is sufficiently low so that high heat generation and high wear are prevented. 
However, contact or close clearance must be maintained at all operating conditions. 

Seal balance can be achieved, at least theoretically, by properly adjusting the secondary 
seal diameter (fig. 8). A common term used by seal designers is the geometric balance 
ratio (or modulus). This modulus is defined as the ratio of the hydrostatic closing area 
Ars to primary seal-face (dam) area Agjj and is used to determine the location of 
the secondary seal diameter. It is desirable to predict this location analytically. 

Unfortunately, a fluid-film seal may be "balanced” only at one combination of 
operating conditions. Balancing is strongly dependent on the variation in film thickness 
and the pressure profile load factor. For gases the pressure profile factor varies with 
the sealed-gas pressure differential. The pressure profile factor is defined as the ratio 
of the net or average seal-face pressure to the sealed-fluid pressure differential. Here- 
inafter, the pressure profile load factor will be referred to as the load factor. Both the 
load factor and the geometric balance ratio have other names in the literature (refs. 1 
and 6 to 9) and are sometimes defined in slightly different ways. Some of these names 
are listed in table I. Since the load factor can equal the geometric balance ratio at only 
one set of operating conditions, it is impossible to completely balance an ordinary face 
seal for all situations. Engineering judgment must be employed to select the proper 
design. 

The importance of the load factor and the geometric balance ratio can be illustrated 
by considering a face-seal force balance. The basic equation defining seal hydrostatic 
closing force is (fig. 8). 
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Net hydrostatic closing force = Fg ± (Fp + Fj) + Ajjg AP - W ^ P dx 

Design philosophies differ; however, a common pressure-balancing practice for gas 
film seals is to select the spring force Fg to overcome only the frictional forces Fp 
and the inertial forces Fj. (The frictional forces are due to the secondary seals (e. g. , 
O- rings and piston rings) and the antirotation lugs (e. g. , torque pins) rubbing on the 
housii^. ) 

A fundamental consideration in designing pressure-balanced seals is the selection 
of the secondary seal diameter. This diameter determines the hydrostatic closing 
force, as illustrated in figure 7. By proper positioning of the secondary seal diameter, 
the closing force can be equal to the seal opening force or, at least theoretically, to any 
degree of seal-face loading. The secondary seal diameter can be found from the geo- 
metric balance ratio, where 


Geometric balance ratio = 


-^HS 

^SD 


Another important parameter is the pressure profile load factor F, which is de- 
fined as the pressure (pneumatic) opening force normalized to the sealed-pressure dif- 
ferential force acting over the entire seal-face (sealing dam) area or 


F 


Pressure opening force 


AP A, 


SD 


Across a seal face. 


Opening force = W P dx 


r 

Jo 


If the seal opening force is equated to the hydrostatic closing force 




W/o Pdx^APAjjg 


Substituting this condition into the load factor relation results in 


— ^HS 

F = —2= = Geometric balance ratio 
■^SD 
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When this situation exists, that is, when the load factor is equal to the geometric balance 
ratio, the seal is said to be perfectly balanced . 

Once the load factor is known, the seal balance diameter can be simply calculated. 

Seal balance diameter = F(2R2 - + 2Rj^ 

(For a perfectly balanced seal, the seal balance diameter equals the secondary seal di- 
ameter. ) For some cases the seal opening force can be evaluated analytically, and 
hence the load factor can be predicted analytically. 


LAMINAR FLOW 


Flow Between Parallel Plates in Relative Motion 


The case of parallel flow (fig. 6) of an incompressible viscous fluid is probably the 
simplest exact solution of the Navier-Stokes equations of motion. However, it is a 
reasonable approximation to many sealing flows. Since the flow is incompressible and 
assumed to be isothermal, p and p will be constant; also the body force is neglected. 
Since the flow is only in the x-direction, solutions are sought where v = w = 0. Sidewall 
effects or side leakage effects are neglected. From the continuity equation we find 
3u/9x = 0, so u = u(z). Thus, the governing equation for this flow is 


dx 


= M 



(20) 


Physically, this means that the pressure drop is proportional to the viscous friction, as 
expected. This equation is integrated twice with respect to z, yielding 



lax 2 ^ 



(21) 


Now apply the boundary conditions (which are the "no slip" boundary conditions since 
the fluid wets the solid wall) 


u = 

at 

N 

It 

O 

u = U2 

at 

z = h 
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Successive substitution into equation (21) yields two simultaneous equations for the con- 
stants of integration Cj and C 2 - The result is 


C. (Uy - U.) 

^ h ^ ^ 2 dx 


^2 “ 

These constants are substituted into equation (21), giving 


2/j, dx 


(U2 - ^l) 

h 


z + Uj^ 


( 22 ) 


It is useful to examine the situation when U|^ = 0. Thus, letting Ug = U, 



This is known as the general Couette flow velocity distribution equation. For the case 
U = 0, that is, when both surfaces are at rest, equation (23) gives a parabolic velocity 
profile such as that shown in figure 9. This is known as plane Poiseuille flow. The 
general behavior of the velocity profile is shown in figure 10. For 3P/9x = 0 and 
U # 0, the velocity profile is linear and is known as simple Couette flow. In general, 
the velocity profile is made up of both the pressure flow component and the shear flow 
component. Also depending on the value of dP/dx, the profile may be convex or con- 
cave, or it may even reverse itself. The point at which the flow will occur in a direc- 
tion opposite to U can be calculated by letting du/dx = 0 at the stationary surface. 

The result is that, for dP/dx > 2U/i/h^, the flow will reverse itself over certain values 
of z. (This type of velocity profile characterizes flow in a viscoseal. ) For incompres- 
sible flow the velocity profile does not vary with axial flow distance. Also since the 
governing equations are linear, the pressure flow and shear flow components can be 
analyzed separately and then superimposed. 

Leakage flow rate . - The continuity equation in differential form can be replaced by 
the condition that the volume of flow in every flow cross section must be constant. 
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- dz = Constant 
h 


(24) 


Q = Wh 



Physically, this means the flow cross-sectional area is multiplied by the average veloc- 
ity. Hence, for this parallel-flow case 


Q = W f udz = W f f— — (z^ - hz) + (Uo - Uj) - + Ujl dz (25) 

After integration and substitution of limits the volumetric flow rate becomes 

Q_:^dP^Wh(„ u.) (26) 

12|U dx 2 ^ 

Note the strong cubic dependence of the leakage flow rate on the gap h. If the gap 
were doubled, the leakage would increase eight times for laminar flow. Although there 
is lesser dependence on gap for turbulent flow and in labyrinth seals, tight control of the 
flow cross-sectional area is essential. The clearance between the seal and shaft that 
forms the flow path should be as small as possible. However, the minimum clearance 
possible is limited by shaft deflections, variations in the film thicknesses, fabrication 
and assembly tolerances, unequal expansions, and so forth. Because of variations 
present in any seal, the gap in equation (26) can be the "mean effective separation" be- 
tween the sealing surfaces. 

Pressure distribution . - The pressure distribution is still unknown but now can be 
found. Rearrangii^ terms and integrating with respect to x yields 

X 

dx + (U2 + Uj) / dx 

Wh^ Jq 

= _ l^ + Ji^(U„ + Ui) (27) 

Q 2 “ ^ 

Wh"^ Wh^ 

Therefore, at any point x in the sealing passage 

p = Pj + 2^ + U2 + Uj^ (28) 

h^W ^ ^ 



le 



at the end of the passage (x - L, P = P2). Hence, the volume flow rate can now be found 
in terms of the usual known physical quantities in fluid-film sealing. Thus, 

3 

^ " TT- ^^1 ■ ^2^ ■" - (^2 + Ui)W (29) 

12/iL ^ 2 ^ 

Comparing equation (29) with (26) shows that the constant pressure distribution is 
merely 


dx 



(30) 


for plane Poiseuille flow. Now the static pressure at any point in the flow is found by 
substituting equation (29) into (28). The result is 


P=Pl-i(Pl-P2) 


(31) 


There is no pressure distribution associated with shear flow. Also the pressure is in- 
dependent of the gap and the fluid properties. 

Opening force. - The seal opening force is found by integrating the pressure over 
the surface area on which it acts: 


F = 




(32) 


F = PjLW - (P^ - ?2) (33) 

2 

This expression defines the force represented by area abed on the P-x diagram shown 
in figure 11. In seal work it is generally of interest to know the increment of lifting 
force due to the pressure differential rather than the absolute pressure force. This 
force is represented by area dec on the P-x diagram and is equal to 


Fnet = F - P2LW 


(34) 


or 
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( 35 ) 


This expression indicates that the effective pressure tending to separate the flow bound- 
aries is one-half of the full pressure differential (linear pressure drop). 

Frictional horsepower requirements . - To determine the power dissipation or fric- 
tional horsepower requirement, recall that the assumption of a Newtonian fluid again 
implies 


9u 
IJ- — 
dz 


(36) 


We will consider two cases. The first case is for a ring sliding on a rod or in the bore 
of a cylinder. The velocity gradient term in equation (36) is found by taking the partial 
derivative of equation (23) with respect to z; 

ini (37 

3zjz_h 2p dx h 

For convenience of manipulation, Uj is taken as zero. The viscous drag force resist- 
ing the sliding motion between the two boundary surfaces is equal to the product of the 
shear intensity and the area over which this acts. Frictional horsepower is in turn the 
product of this drag force and the relative slidir^ velocity (for a ring or bore configura- 
tion, W = ttD). Hence, 


?tDUq 

HP = L T dx 

6600 


(38) 


Combining equations (30), (36), and (37) and substituting in equation (38) yield, after in- 
tegration, 


6600 ^ \2 L 



dx 


jrpDLuf 

HP = = 

6600h 


ffDhUj 

13 200 


(Pi 


P9) 


(39) 


(40) 


Another important case is a clearance-type rii^ seal where the boundaries are ro- 
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tating with a relative angular velocity Q in the y-direction (Uj = U 2 = 0). Only the ve- 
locity gradient in the y-direction applies. Here we may assume a linear velocity gradi- 
ent in that the fluid at each boundary has zero velocity relative to that boundary (no 
slip). This is a simple Couette (shear) flow profile. 



(41) 

2 h 

3y 2h 

(42) 


For this case the frictional horsepower becomes, after evaluation, 

HP = 

26 400 h 


(43) 


Center of pressure . - The center of pressure in the flow passage direction can be 
found from 


X„ = 


/ 


P - P. 


W 


X dx 


(44) 


^c = 




(45) 


^c=- 
^ 3 


(46) 


By examining the triangular pressure gradient in figure 11, this result is expected. 
This is the value of the centroid for a triangle. 

Axial film stiffness. - The axial film stiffness is found from 


K = - 


dz 


K s 0 for a parallel surface 


(47) 

(48) 
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Since the seal opening force does not depend on film thickness for parallel films, there 
is no axial film stiffness. Some auxiliary device, such as self-acting lift pads or a ra- 
dial step, must be used to obtain axial film stiffness. 

The procedures for converging and diverging surfaces and for radial and cylindrical 
flow are identical to the procedure for the parallel-surfaces case presented. However, 
the mathematics is more complex because of the geometry variation. 


Flow Between Flat Converging or Diverging Plates 


The governing momentum equation 


^ ^ 


0X 




(49) 


and velocity distribution 


2(1 dx h 


(50) 


are the same as for parallel surfaces. The film thickness variation for converging or 
diverging surfaces can be expressed as (fig. 12) 


h = hj + ax 


Leakage flow rate . - The leakage flow rate can be found from 


h?W 


Whn 


Q = - _L_ if (1 + ax)-^ + U«(l + ax) 

12(1 dx 2 ^ 


(51) 


Rearranging terms and integrating with respect to x yields 


P = P 


6 (ih 


hi(h2 - hj^) 


_S. 


2 .2 \ ^ 1 ^ 


Lhih^ 




(52) 


Thus, the volumetric leakage flow rate per imit length is 
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(53) 


Q (Pi - P2> , '■l(^ 

W 6/iL (^ + 1) (/3 + 1) 


(U2 + Uj) 


where 



When /3 = 1, we have parallel plates 


h^W 

12juL 


(Pi - P9) + 


h^(U2 + U^)W 


Note that 


/3 < 1 for converging plates 
|3 > 1 for diverging plates 


(54) 


Generally, tapers in seals are small, but in some applications it is desirable to 
machine a converging taper in order to achieve a positive film stiffness. Also it can 
readily be shown that 


, 3.2 . 2.2 

^ + 1 2h^ 


(55) 


By defining a characteristic film thickness 


char 


r u2i,2 

h^ 

^m 


.1/3 


(56) 


Thus, the leakage flow rate can be expressed as 


Q = 


■■Lar^<Pl - 

12 pL 


+ 


Whj/3(U2 + Uj) 

(/3 + 1 ) 


(57) 
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When there are no tilts present, hj = h 2 (hence, a = 0). Then = h and the 

classical parallel- surface mass leakage equation is evolved. However, if an effective 
or linear apparent tilt is desired, the effect on mass leakage can be easily calculated 
from the parallel- surface leakage equation. A simple computation enables leakage with 
deformation to be readily calculated or read from a parallel-surface leakage plot. 

Pressure distribution . - For the case where the relative surface velocity in the flow 
direction is zero, the pressure distribution across the sealing dam is 

P = Pj - (Pj - P 2 ) - ^ (58) 

L 


where the linear tilt factor ^ is defined as 

h?(2hi + ax) 

3' = — (59) 

2hj^(hi + axf 

For parallel surfaces, a = 0; hence, = hj = h 2 and the tilt factor ^ is 1, which 
yields the parallel- surface pressure distribution. Representative pressure profiles for 
parallel, converging and diverging surfaces are shown in figure 13. 

Opening force . - Upon integration of the pressure distribution, the net opening force 
can be found and is 


^ (=>1 - => 2 ) ^ ^ « l )[“ - 

1 + /3 (g _ L /3 + 1 J 

^ ^ ' Y ' 

Hydrostatic Hydrodynamic force 

force 

When there is no relative motion of the surfaces in the leakage flow direction, the net 
surface opening force is 


F 


net ■“ 


WL AP 

1 + P 


The load factor is 


'Ago AP" 1 +^ 
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For parallel surfaces, /3 = 1 and 


F = 0. 50 

For converging faces, /3 < 1; and, for example, if j3 = 1/2, 

F = 0.67 

For diverging faces, /3 > 1; and, for example, if /3 = 2, 

F = 0. 33 

The practical range of load factors appears to vary from 0. 6 to 0. 9 for liquid face seals 
currently in use. 

For /3 = 1, the parallel-surface case, equation (60) reduces to 

^■net = ^ - Pj) (35) 

The second term is the load-carrying capacity of a slider bearing which is infinite in 
y-direction. 

When = ? 2 , the force induced in the classical hydrodynamic slider bearing Uj 
is equal to zero. This pressure distribution is shown in figure 14. The dimensionless 
hydrodynamic force is defined as 


^net 

PaWL 




Pah 


2 

1 


where A is the bearing number. The bearing number is similar to the load factor ex- 
cept that the ambient pressure is the reference pressure. Hydrodynamic pressure gen- 
eration is discussed in the section MECHANISM OF FILM PRESSURE GENERATION 
BETWEEN RELATIVELY MOVING SURFACES, 

Frictional horsepower requirement. - First, we consider friction on a moving sur- 
face, for example, a ring sliding on a rod or in the bore of a cylinder. From 

2 ax h 


we get 
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HP = 


iraDhjU2(Pi - P2) 


(61) 


rriJ.. 


DLU^ 


3300(^ - l)hj 


^2 ln/3 - 


3^^ 
p + V 


6600()3 + 1) 


Second, for clearance-type seal rings, boundaries are rotating with relative angular ve 
locity S2, and = U2 = 0. Only the velocity gradient in the y-direction applies; that 
is. 


w = 



1 

h 


The result is 


HP = J£J_ 

26 400hj /3 - 1 


(62) 


Radial Flow Between Annular Disks 

The governing equation for the radial flow between annular disks, representative of 
mechanical face seals (fig. 15), is 



3P a^u 

- — + /I - — 


0r 


3z" 


(63) 


The radial geometry yields a centripetal inertial force (centrifugal force) because of the 
curvature on the left side of equation (63). This equation yields the following velocity 
profile: 


u = 


(z^ 

12/ih^ 



2/i 9r 


- hz) 


(64) 


Leakage flow rate . - Use the continuity equation for volumetric flow in any radial 
direction r 


.. ^ pn^Wrh^ 

Q.wVo udz=^^ 


wh^ ap 

12 p ar 


(65) 
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Let 


W = 2vr 


Integrating with respect to r between Rj and r yields 


P = 




At r = R 2 , P = P 2 ; hence, the flow rate can be found from 


Q = 





( 66 ) 


(67) 


Q = 0 


when 


r2 . r2 , -20^ (P 2 . pp 

2por 

When AP > 0, leakage increases due to centripetal force; when AP < 0, leakage de- 
creases due to centripetal force. 

Frictional horsepower requirements . - The frictional horsepower can be found by 
using 


r = = 

dz h 


Thus, 


HP = r^2 J.3 

6600h -^Rj 


(69) 


or 
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( 70 ) 


HP = _£i!«i. (r4 
13 200h ' 



Opening force . - The opening force can be found by integrating the pressure distri- 
bution equation 


P = 




(71) 


where 



For Pj > ^ 2 ’ ^net *^han for inward flow (when Pj^ < P 2 )- For the hydrosta- 

tic case and R 2 /Rj^ !• 1> it can be shown that the static pressure is linear and 

F„et - f (Pi - P2 )(r 1 * "l'*2 - 2K?) 


Pressure distribution . - The pressure distribution due to the centrifugal force can 
be examined in more detail by solving for pressure distribution when the inner and outer 
diameters do not have an imposed pressure differential (Pj^ = P 2 ). For this case, equa- 
tion (71) becomes 


P = 




(75) 
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For most face seals, R 2 /R 1 is close to 1. Hence, 



Ro - Ri 


Ri 


In 


r 


R 


1 



(This is known as the face-seal approximation. ) 


P = Pi - P^hr - R 1 KR 2 - r) (76) 

Equation (76) indicates that the centrifugal force term is always negative and inde- 
pendent of the nature of pressurization. For both internal and external pressurization 
the centrifugal effect diminishes the pressures on the seal face, reducing the separation 
forces and resultii^ in a decreased seal clearance. 


Axial Flow Between Annular Cylinders 

Concentric cylinders . - The governing equation for the axial flow between concen- 
tric annular cylinders, representative of bushing and circumferential seals (fig. 16(a)), 
is 


^ + 1 


(77) 


Note the additional friction term due to curvature. The velocity distribution can be 
foxmd from 



Upon integrating u between and R 2 , the flow becomes 
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( 79 ) 



For the special case of a circular cylinder (Rj^ = 0), the velocity distribution is 
found from 


u = 


4/x dz ' ' 



The maximum velocity is 


^max 


1 dPp2 

4p dz ^ 


Since the velocity profile is parabolic, 


u 


av 


^max _ ^ jj2 

2 8/i dz ^ 


(80) 


(81) 


(82) 


Hence, 




2 


8pL 


(83) 


Eccentric cylinders . - Eccentric cylinders simulate the case of two nonparallel de- 
veloped surfaces. Let h = hj^^(l + c cos d). Thus, 


u ^ y[y - hjj^(l + e cos 0)] = f(y, 9) 

2/uL 


(84) 


Now 


n h p 


y[y - h^(l + € cos 0)]dy R d9 


(85) 
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(86) 


;rPiRh;L / 
Q = _i. ^ (i 

6/liL V 



The eccentricity can be due to loading, fabrication, and/or assembly tolerances. 

The actual eccentricity of the shaft and bushing is very seldom known; therefore, only 
the two limits, maximum and minimum flow, are usually computed. For e = 0, equa- 
tion (86) reduces to equation (79). Comparing these two equations shows that the flow 
rate will be 2^ times that for a concentric cylinder when the cylinders are fully eccentric 
(e = 1). An exact analytical treatment is complex. Details and a partial solution can be 
found in reference 10, 

The higher leakage for the eccentric cylinders can be thought of as being due to a 
larger effective gap and, hence, a larger shear surface area. 


TURBULENT FLOW 

For turbulent flow, exact physical knowledge is unknown. Hence, an exact differ- 
ential analysis model such as the one that describes laminar flow is impossible to solve 
or impractical for design analysis purposes. Thus, approximate solutions must be 
found. A widely used method in fluid mechanics and hydraulics is the approximate inte- 
grated average method. Although the integral models only satisfy mean conditions in the 
flow field, they have shown good results on gross quantities such as seal leakage and 
pressure distribution. However, it will be seen that an empiricism is required to find 
a solution. 

The viscous friction is balanced by the pressure drop. This is the classical fluid- 
flow case. This model is widely used to describe pipe and duct flows. 

Consider the control volume shown in figure 17 for situations when the fluid inertia 
is negligible. The momentum conservation is a balance between the pressure and the 
viscous friction force, which is 


A dP = -T^ dA^ 


(87) 


Introducing the parameters hydraulic diameter D 

D = 

dx 


and mean Fanning friction factor f 
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into equation (87) results in 


f = 


T 


W 



2 


D— = -2pu^f 
dx 


(88) 


Substitutii^ the mass flow definition 


M = puA 


(89) 


yields the following useful form: 


dP = 


-2M^ dx 
pDA^ 


(90) 


The pressure at a position x along the leakage length can be found by integrating 
from the inlet to any position x 


f-f 


2M^f 

pDA^ 


dx 


(91) 


or 


P= 


2M^fx 

pDA^ 


(92) 


At the seal exit, x = L and P = P 2 ; hence, the mass flow rate can be found from 


M = 


fp DA^(P^ - Pg) 
2fL 


(93) 


Substitutir^ equation (93) into equation (92) results in 
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(94) 


P=Pl-(Pl-P2>f 

Hence, for either laminar or turbulent flow, the pressure distribution is independent 
of the fluid properties and the film thickness. For radial flow between coaxial parallel 
disks and parallel plates, the hydraulic diameter D is given by 

D = 2h (95) 

Generally, the mean friction factor is related to Reynolds number by a expression of 
the following form: 


(Re)" 

It is useful to express the Reynolds number as 


Re = 


2M 

Wp 


(96) 


(97) 


For laminar flow the friction factor is derived from the classical, viscous flow solution 
(eq. (23)) and the derivation is also shown in reference 11. The resulting mean friction 
factor - Reynolds number relation is 


24 _ 12jiW 
Re M 


(98) 


Thus, 


. phV(P. - P„) 

M = 1 £_ 

12pL 


(99) 


This is the identical result obtained in equation (57) by setting Uj = U 2 = 0 and M = pQ. 


Turbulent Flow Between Parallel Plates 

The Blasius relation of friction factor and Reynolds number appears to satisfactorily 
describe a large class of fully developed flows even though it is experimentally deter - 
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mined. Thus, in equation (96), k = 0. 079 and n = 0. 25. Substitution in equation (93) 
yields 


21/7 4/7,12/7^( )4/7 

M = 1 -1 2 (100) 

(0.079)^/'^L‘^/'^fi^/'^ 

which gives the functional relation of the variables in quasi-fully-developed turbulent 
flow. The leakage dependence on film thickness is no longer cubic, but it is less than 
quadratic. However, the leakage is still most sensitive to gap. The fimctional variation 
of leakage is compared for laminar and turbulent flow in table n. 


Turbulent Flow with Small Linear Deformation (and Constant Width) 

Again, as in laminar flow, the surface deformation is represented by 

h = hj^ + ax (101) 

The mass leakage flow rate is identical to the parallel surface equation (eq. (100)), ex- 
cept that the characteristic film thickness is used: 

M (102) 

( 0 . 

The pressure distribution is also identical to that found for laminar flow (eq. (94)). 

The dependence on flow-width is the same as for laminar flow; but there is lesser 
dependence (4/7) on density, flow length, and pressure differential. Also notice the weak 
dependence on molecular viscosity - a characteristic of turbulent flow. However, tur- 
bulent flow is characterized by large-scale momentum exchange (eddies). This macro- 
scopic fluid behavior can be represented as apparent shear stresses. Hence, a turbulent 
viscosity can be orders of magnitude higher than the molecular viscosity. Thus, turbu- 
lent flow leakage is characterized by a high effective viscosity, which means lower leak- 
age than predicted by laminar-flow models but also higher shear heating. 


Turbulent Radial Flow 

Equation (90) can be integrated for radial annular flow without relative seal-face 
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rotation by using the same restrictions as previously stated. The following result is ob- 
tained for mass flow (here W = 2nT and x = r): 


M =. 


2.p4Ahl2/7(p^ _ p^)4/7 

/ \4/7 

(0.079)^A^1A/_J L 


(103) 


The pressure distribution can be found from 


(Pi - P2)/ 


P = Pl- 


(R 


3/4 3/4 


(104) 


Rf4 r3/4 

The leakage flow rate with rotationally induced turbulence (p. 61) is 


M = 


3gp 


1/4 


1 


4(0-0195)\^p3;,fi3y (R3/4.R3/4j 


jPj - P 2 + (0. 128)p^R|f2^ - 


(105) 


Sneck (ref. 12) used an eddy viscosity approach to derive equation (105). Details 
can be found in that paper. 


Turbulent Flow Between Concentric and Eccentric Cylinders 

When flow is turbulent, leakage is not affected as much by eccentricity as when 
flow is laminar. Generally, turbulent flows are less sensitive to gap, as shown in ta- 
ble n. This is probably due to the large momentum exchange that characterizes turbu- 
lent flow. Leakage for the fully eccentric shaft is 30 percent higher than for the per- 
fectly concentric shaft. The analytical expression for turbulent flow through an eccentric 
annular gap is somewhat complicated and approximate and is given in graphical form in 
figure 18 (ref. 10). Figure 19 compares the effect of eccentricity on laminar and turbu- 
lent flow. 
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APPLICATION OF THEORY TO FACE -SEAL CONCEPTS 


An application of the preceding derivations will now be cited. References 13 and 14 
numerically solved six face-seal-contour cases on a digital computer, relating film 
thickness, leakage flow rate, opening force, and exit pressure. These contours repre- 
sented converging, diverging, and symmetric faces and were compared with parallel 
faces. The studies of references 13 and 14 were conducted with sea water at 300 K 
(80° F) as the fluid medium. The seal inside diameter was 22. 162 cm (8. 725 in. ) and 
the outside diameter was 24. 067 cm (9. 475 in. ). The face contours are shown in fig- 
ure 20 (from ref. 14). 

The relation of minimum film thickness to inlet water pressure is shown in figure 21 
for a leakage flow rate of 0. 002 m^/min (0. 5 gal/min) and a change in film thickness Ah 
of 0. 00127 cm (0. 0005 in. ). All nonparallel cases resulted in a decrease in minimum 
film thickness from the parallel case. (The centrifugal force had a negligible effect on 
the film thickness. ) 

Figure 22 shows log-lc^ plots of the minimum film thickness against the ra- 

tio of leakage flow rate to inlet water pressure, resulting in a single curve for each face 
contour. 

The radial pressure profiles illustrated in figure 23 show the effect of the various 
contours. A rapid dropoff in pressure occurred near the outer radius with the diverging 
surfaces; a slow decrease occurred with the converging surfaces. 

The resulting plot of seal-face load against inlet water pressure is shown in fig- 
ure 24 for a constant leakage flow rate and the same inlet and outlet film thickness dif- 
ferences. Only the symmetrical parabolic contour is similar to the parallel case; the 
others differ by an increasing amount as the pressure increases. The effects of Ah 
and Q on the seal-face load are shown in figures 25 and 26, respectively, for a con- 
stant P 2 * Only the symmetrical parabolic face load remains constant as both Ah and 
Q are varied. All approach the load of the parallel case as Ah approaches zero with 
constant Q. However, the diverging case approaches zero and the converging case ap- 
proaches the full load of 1. 43x10^ N (32 150 Ibf) as Q approaches zero with constant 

Ah. 

Finally, figure 27 illustrates pressure profiles for a converging gap in the flow di- 
rection, for parallel surfaces, and for a diverging gap, including the effect of a pro- 
nounced area change. 


COMPRESSIBLE FLOW OF GASES 

The equations previously derived have been for liquids, which can be considered to 
be incompressible under most sealing conditions. The constant-density assumption 


34 



means that the equations will be simpler to solve and that pressure appears only in a 
dynamic role not thermodynamically. Hence, the previous solutions were obtained with- 
out use of the energy equation. This is not the case for compressible flow. Since den- 
sity is now a dependent variable, the momentum equation is interrelated to the energy 
equation. Also another equation is necessary, the equation of state. The physics of the 
flow are also different and more complex. Principles of gas-film-seal fluid mechanics 
are discussed in great detail in reference 11. Herein only the salient features will be 
pointed out. 

Generally, for thin films the standard lubrication assumption of isothermal surfaces 
can be made even for gases. This is so because gas lubricant films are very thin and in 
contact with relatively large masses of metal. However, at higher flow rates, an adi- 
abatic or polytropic equation of state may be more relevajit. This isothermal or poly- 
tropic equation of state enables a relation to be found which eliminates consideration of 
temperature in solvii^ the set of equations. Hence, the energy equation need not be 
considered. The polytropic equation of state is 

p 

— = Constant 

The gas film theory, as well as gas dynamics, is based on the premise that at any point 
in time and space the flow is in thermodynamic equilibrium and the perfect-gas law ap- 
plies 


P = psn 

In summary, it is the local change in density in the film that characterizes com- 
pressible flow. Because of the isothermal assumption, the density change results pri- 
marily from changes in pressure as predicted by the perfect-gas law. Hence, we expect 
to see differences from incompressible flow in physical quantities of interest. Also the 
viscosities of gases increase slowly with temperature, as opposed to the rather rapid 
decreases that liquids exhibit. This is discussed in more detail in appendix B. 


Compressible Flow Equations 

The set of equations that will be used for compressible fluids (gases) between nar- 
rowly spaced annular disks in polar coordinates (fig. 28) is 
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Continuity: 


(106) 


- — (pru) + — (pv) = 0 
r 0r dz 


Radial momentum conservation; 




3r 


dz 


dr 


dz^ 


(107) 


Circumferential momentum conservation; 


a2 

9 w 


= 0 


9z^ 


(108) 


Axial momentum conservation: 


^ 0 (109) 

0Z 


Equation of state; 


P = p^PT 


( 110 ) 


Energy: 


pcuil.pcvil = ki!l 
P 3r P 3z 5^2 


( 111 ) 


The formulation of the appropriate equations governing gas-film -seal behavior is 
detailed in reference 11. 


Classical Viscous Compressible Flow Model 

Parallel radial flow . - The classical hydrostatic radial leakage-flow formula illus- 
trates the effects of compressibility. It is valid for compressible viscous flows and is 
the widely used gas leakage formula in the seal industry. M rotation effects are ne- 
glected and the flow is assumed to be isothermal, the governii^ equations reduce to the 
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equations for mass and radial momentum conservation. The radial momentum equation 
is 


dP _ afu 


dr 


3z" 


( 112 ) 


This yields 


u = — ^ + C.z + C„ 

2/i dr ^ ^ 


(113) 


Upon applying the boundary conditions the result is 


u=-L®(z2-hz) 

2/i dr 


(114) 


The mean velocity at any radial cross section is 


_ 1 
V - 7 


dz 


(115) 


v = - 


dP h^ 


dr 12p 


(116) 


Now, the mass leakage flow rate is 


M = pWhu 


m 


= /_ pWh^ VdPX 

\ 12p/drj 


(117) 


If the perfect-gas relation is substituted and equation (117) is rearranged, the radial 
pressure gradient can be expressed in the following forms; 


dP ^ _ 12p^TM 


dr 


(118) 
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p dP = - jjj. 

Wh^ 


(119) 


Since 


W = 2ht 


jjp2 ^ 12 m ^TM dr 

' ' .h3 - 


or 


p2 = - 12^ ^TM In r Q 


Applying the boundary condition at r = Rj, 



( 120 ) 


yields 

p2 = P?+i2ii^ln^ 

^ _u3 r 


Applying the boundary condition at r = R 2 » 
yields the classical mass leakage flow equation 

^1 

12p^?T In-Jt 

«2 


( 121 ) 


( 122 ) 


Using equation (122) in equation (121) yields the following pressure distribution equation: 
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11/2 


P=Pl 1- 


,2\ln-i 


p2 1 Ri 


The classical solution for the pressure distribution is independent of both the fluid prop- 
erties and the gap. Also, the velocity profile does vary with radius, and hence for com- 
pressible fluids the flow is not fully developed but is quasi-fully developed. 

Face-seal approximation (AR/Rj « 1). - Ih many gas film seal faces, AR is much 
less than Rj. For this case, ~ 


In^^-M 


By setting W = 2irr, equation (122) becomes 


Wh3(p^ - 4) 

24]Lt^T(R2 - Rj) 


Equation (124) can be written as 


Since P^^ = (Pj + P2)/2, 


Wh'^(Pj + P2>AP 
24|Li^T(R2 " 


L = R2 “ 
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M = 


Wh^p^y AP 
12pL 


The preceding equation is written in incompressible form. 

Leakage expressed as standard cubic meters per second: In order to have a mean- 
ingful assessment of leakage at a common base (standard conditions) to account for den- 
sity variation with temperature, it is convenient to express the leakage in the following 
manner: 


Q 


scfm 


M 

^std 


^std 


Therefore, 


Q 


scfm 


Pstd , 



and equation (123) becomes 


P = 


1 - 



(125) 


where 

X = r - 

(Mass flow is proportional to Pj - I >2 for compressible flow but proportional to the 
pressure drop for incompressible flow. Also compare eq. (125) with its incompressible 
counterpart, eq. (31). ) This is the classical form of the narrow-slot leakage equation; 
thus, hydrostatic radial flow can be approximated with a narrow-slot (plane flow) anal- 
ysis with small error if (R 2 ~ Rj^) This approximation shall be known, herein- 

after, as the "face- seal approximation. " 

Opening force: The net seal opening force is found by integrating equation (125). 

The result is 
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(126) 


L 

w 



By algebraic manipulation, it can be shown that the load factor 



This equation is plotted in figure 29 and illustrates the strong variation of load factor 
with pressure ratio. Since F = F(P 2 /Pj), F can equal the geometric balajice ratio only 
at one operating condition. This illustrates why it is impossible to completely balance a 
face seal for all sealed-gas pressure situations and that a slight unbalance at large pres- 
sure differential will greatly overload the seal. 

Since equation (126) also describes the pressure profile for laminar and fully de- 
velqjed turbulent flows, the preceding load factor applies to turbulent flow as well. (The 
load factor equation is independent of film thickness and fluid properties. ) 

For pressure ratios close to 1 the gas should behave like an incompressible fluid. 
This is the case, since the load factor is 0. 5, At very large pressure ratios, 


2 

3 


However, this high-pressure-limiting situation is not physically correct. In reality, the 
governing equations break down, in that fluid inertial forces become important and the 
phenomenon of fluid choking can occur. This is discussed in the section Physics of com- 
pressible flow and concept of choking . 

Center of pressure: The center of pressure can also be readily foimd from equa- 
tion (126), which yields 
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(127) 


W(Ro - 



.A 

Vj 15 



F 


\pii 

)] 

2 


min 


> 


An alternate formulation of this classical compressible viscous flow case would ave 
been to start with the compressible Reynolds lubrication equation (ref. 15). 

nn^Ri-fully-develoned flow. - Incompressible flow in long narrow spaces is c arac- 
teriz^byalinearpre1i;i;r^op and a velocity profile that does not change in the low 
direction. This can be seen by examining the continuity equation for this parallel flow. 


3 ^ = 0 

9x 


Let us now examine the continuity equation for the same case but for a compressible 
fluid: 


9(pu) = 0 

0X 


or 


Iksi = + u.^ = 0 

9x 0x 0x 


Hence, from the equation of state, P « p; that is, as the pressure drops the density also 
decreases, and the velocily Is no longer constant in the flow direction but mcreases 
(tig. 30). Thus, the term ”quast-fuUy-developed flow" is used in compressible viscous 

flow 

A gas generally has a large specific volume increase associated with a pressure 
drop. (Liquids are, of course, incompressible. ) Hence, for a given pressure op 
across a seal, the less dense the sealed fluid, the lower the resultant mass flow rate^ 
(Perhaps this is one of the reasons that a labyrinth seal is effective in sealing a gase 
fluid but would be ineffective for the same fluid in the liquid state. ) This lower mass 
flow rate can be seen by comparing the leakage rates for the same gas and through 
same geometry but obtained at two different temperatures (fig. 31, from ref. 16). 
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Pressure Distribution - Application Example 


Figure 32 shows, for parallel sealing faces, the normalized pressure as a function 
of radial position within the sealing gap. The normalized pressure distributions for 
three operating design points representative of aircraft climb, cruise, and idle are 
shown for both compressible and incompressible fluids. Since the faces are parallel, 
the profile is independent of film thickness. However, the profiles shown are valid only 
below the choking points (Mach number, <0. 845). The compressible pressure profile is 
parabolic rather than the linear profile as is found for incompressible fluids (and for 
flow between parallel surfaces). 

Nonparallel surfaces can have been deliberately machined as a taper or can be 
caused by distortions. Distortions of the primary seal faces are inherently present in 
gas film face seals. Distortions present include radial and axial displacements due to 
the centrifugal force and are especially important under high rotational speeds as antic- 
ipated for advanced aircraft operation. A typical centrifugal distortion is shown in fig- 
ure 33. Another common face distortion is thermal coning, which is caused by an 
axial thermal gradient along the shaft. The hotter end of the shaft causes a differential 
shaft radial displacement which results in the face coning illustrated in figure 34. aher 
distortions could be caused by (1) pressure - due to high pressure drops and improper 
seal balance diameter; (2) mechanical effects; (3) asymmetry of rotating seal seat; and 
(4) tolerance buildup due to fabrication and assembly. Generally, for internally pres- 
surized seals, the distortions will cause seal faces to diverge. 

For nonparallel sealing surfaces, the pressure gradients are dependent on film 
thickness. ^ This is shown in figure 35 for the design point of 45 N/cm^ abs (65 psia) and 
311 K (100° F). The pressure gradients are plotted for mean film thicknesses of 5. 1 and 
10. 2 pm (0. 2 and 0. 4 mil) and for a seal face deformation of ±0. 001 radian. (A negative 
(-) deformation indicates a converging film in the leakage-flow direction, and a positive 
(+) deformation indicates a diverging film. ) The Important point (fig. 35) is that for the 
convergent film, a decrease in film thickness results in greater area (larger opening 
force) under the pressure gradient curve; thus, the convergent film has a positive film 
stiffness. Conversely, as shown by figure 35, a divergent film has a negative film stiff- 
ness. This means that for a divergent film a decrease in film thickness (which always 
occurs dynamically) will cause a decrease in opening force. 


Approximate Models of Compressible Flow 

As was the case with liquids, approximate models can be formulated for gases to 
predict seal behavior where an exact model is impossible to solve or impractical for 
design analysis purposes. Such cases are turbulent flow, where exact physical knowl- 
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edge is unknown, and near-choked and choked flow, where nonlinear behavior charac- 
terizes the flow. Approximate integrated average methods can again be used. 

First, quasi -fully- developed compressible flow models will be presented for both 
parallel and small linear deformed surfaces. Then a constant-area case will be formu- 
lated, which is especially useful for choked flow but also includes the quasi-fully- 
developed flow cases. A variable-area analysis is formulated for radial area expansion 
and small surface deformation in reference 17. This analysis includes both the quasi- 
fully-developed flow and constant-area cases, but the solutions are of the numerical 
Runge-Kutta type and will not be presented herein. 

Quasi-fully-developed flow models . - Substitution of the perfect-gas law and mass 
flow definition into the viscous momentum equation (87) yields the following fundamental 
form; 


p dp - -2f^TM (128) 

Da2 

The leakage flow rate results for various special cases are shown in table HI. The func- 
tional relations among the variables and the mass flow rate are similar to those for in- 
compressible turbulent flow. 

Con stant-area case - critical and supercritical choked flow . - When the radial flow 
is close to choked or is choked at the exit, the viscous flow analysis is no longer valid. 
Fluid inertia, neglected in that analysis, becomes important. The rigorous inclusion of 
inertial terms necessitates the solution of nonlinear partial differential equations even 
for parallel surfaces and is therefore not practical for engineering calculations; hence, 

an approximate analytical model must be used. 

It is assumed that the flow in the seal leakage-flow region behaves as a constant- 
area adiabatic flow with friction. A quasi-one-dimensional approximation is made 
wherein it is assumed that the flow properties can be described in terms of their cross- 
sectional averages. 

The following assumptions are made in the analysis: 

(1) The area expansion due to radius increase is neglected. (In many mainshaft face 
seals, the ratio of inside diameter to outside diameter is about 0. 98. ) 

(2) The flow is adiabatic. 

(3) No shaft work is done on or by the system. 

(4) No potential energy differences are present, such as those caused by elevation 
differences. 

(5) The fluid behaves as a perfect gas. 

(6) The sealing surfaces are parallel. 

With these assumptions, the flow is commonly known as Fanno line flow (ref. 18). The 
governing equations when area chaises are neglected are as follows. 


44 



Conservation of mass: 


which reduces to 


Conservation of energy: 


M = puA = Constant 


P 2„2 


P u 

e' + — + — = Constant 
P 2 


or 


i + — = Constant 
2 

where i is the specific enthalpy. This can be written as 


di = c dT = - - du^ 
P 2 


This equation can be reduced to 


^ H- (2L-. .. 1 ) 

T 2 


M^ — = 0 


u 


(129) 


(130) 


(131) 


(132) 


(133) 


where M is the Mach number, u/y/y^T. 
Equation of state: 


or 


P = P^T 


(134) 


^ = (135) 

P p T 
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Conservation of momentum: 


(136) 


-A dP - dA^ = M du 

After the hydraulic diameter and Fanning friction factor are introduced, the set of equa- 
tions (130), (133), (135), and (136) can be combined into a single equation in terms of the 
Mach number alone. Details of obtaining this equation are given in reference 11. The 
result is 


4f dx = (1 - M^)dM^ (137) 

D 

This equation can be integrated for subsonic, critical, and supercritical flow conditions. 
The equation is in an implicit form. Solutions are found by using a linear iteration 
scheme. Details and examples are given in reference 11. 

Physics of com pressible flow and concept of chokii^ . - Compressible flow and 
choking occur in the following way: The static pressure drops to overcome flow friction. 
This pressure drop increases the specific volume of the fluid. Since the area change is 
negligible, the mean velocity must increase as the specific volume increases in order to 
maintain the same mass flow rate at each section of the leakage path. As the friction 
causes the velocity to increase, the fluid momentum change also becomes important. In 
order to accelerate the fluid to a greater velocity, a force is required which can be at- 
tained only through an additional pressure drop. This results in a still greater increase 
in specific volume. This process will continue until the end of the leakage path or until 
the fluid reaches the maximum (choking) condition. Choking occurs when the Mach num- 
ber is 1 at the exit (i. e. , exit velocity is sonic). Should the Mach number reach 1 some- 
where along the leakage path interior, it is expected that behavior similar to duct flow 
will occur. For duct flow it can be shown that the flow process will adjust itself, until 
the point at which the Mach number reaches 1 is shifted to the exit of the leakage pass- 
age. The mass flow rate of the fluid is the maximum which can be handled for a given 
inlet density and passage cross-sectional area. This flow process is known classically 
as Fanno line flow. 

The entrance velocity head loss is negligible in subsonic viscous flow compared to 
the total seal pressure drop; however, under choked conditions the entrance velocity 
head loss is no longer negligible. This results in an entrance pressure loss. Also, 
under choked flow conditions, the exit pressure is larger than the sump pressure and in- 
creases with film thickness. The pressure decreases throi^h e 3 q>ansion waves to the 
sump pressure in the outer cavity, (However, this does not significantly affect the axial 
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force balance on a seal nosepiece. ) 

Application examples. - The computer program from reference 9 was used to carry 
out the quasi-one-dimensional flow analysis. By the use of this program, several prob- 
lems have been solved. One such problem analyzed is a case which simulates advanced 
aircraft engine mainshaft seal conditions. The sealed-gas pressure and temperature 
were 45 N/cm^ abs (65 psia) and 311 K (100° F), respectively. The exit ambient pres- 
sure was 10.3 N/cm^ abs (15 psia). The seal- face radial length (seal leakage passage 
length) was 0. 127 cm (50 mils). Figure 36 shows the calculated sealed-gas leakage as a 
function of film thickness, which ranges from 3 to 51 pm (0. 12 to 2 mils). 

This approximate analysis agrees with the classical differential analysis for com- 
pressible viscous subsonic flow, which is used when the Mach number is less than 
1/Vy. As shown in figure 36, for film thicknesses less than 8 pm (0. 3 mil), the Mach 
number is less than 1/v/y = 0. 845 and the leakage dependence on film thickness is the 
classical cubic dependence. The isothermal viscous flow model loses its validity when 
the Mach number exceeds l/^/y (ref. 18). The present approximate model, however, 
is valid for both subsonic and choked flow. Figure 36 shows that for gaps larger than 
8 pm (0. 3 mil) the leakage has a lesser than cubic dependence on film thickness. (Chok- 
ing occurs at a film thickness of 13 pm (0. 52 mil). ) The limiting case of orifice flow, 
in which the flow rate varies linearly with film thickness, would be achieved when the 
film thickness approaches the order of the seal-face radial length of 0. 127 cm (50 mils). 

As indicated in figure 36, the transition from laminar to turbulent flow occurs for 
a film thickness of approximately 33 pm (1. 3 mils). A Reynolds number, based on hy- 
draulic diameter, of 2300 has been chosen to be the critical transition Reynolds number. 
This appears to be a universal critical transition Reynolds number for flows in ducts, 
pipes, and bearings. For laminar flow a mean Fanning friction factor of 24/Re2j^ was 
used. This friction factor is derived from the classical viscous compressible flow so- 
lution. For turbulent flow the Blasius friction factor - Reynolds number relation was 
chosen. The mean friction factor equals 0. 079/Re For the transition flow regime 
the exact nature of the flow, and hence the friction factor, is complex and not fully un- 
derstood. The friction factor used for this flow regime is derived in reference 9 by as- 
suming a smooth transition from the laminar friction factor to the turbulent friction fac- 
tor. The Reynolds number range of 2300 to 3000 is arbitrarily selected as the 
transition flow regime. 

Correlation with experiment. - Experiments were conducted at the Lewis Research 
Center on a radial geometry representative of face seals to check the validity of the 
quasi-one-dimensional model for choked flow. Leakage flow was studied for two coaxial 
rings, 14 cm (5. 50 in. ) in inner diameter, 15. 2 cm (6. 00 in. ) in outer diameter, and 
separated by a fixed parallel gap of 38 pm (1. 5 mils). The reservior pressure of 41. 8 
N/ cm^ abs (60. 6 psia) was held fixed, and the exit pressure varied. Figure 37 shows a 
comparison of this quasi-one-dimensional analysis, the classical differential analysis 
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for compressible viscous subsonic flow, and the radial flow experiment. The experi- 
ment shows that the flow does become choked. The maximum flow rate (choked) is 
0. 0080 kg/sec (0. 0176 Ibm/sec). The present analysis predicts a slightly higher mass 
flow, but agreement is within 19 percent over most of the range. On the other hand, 
classical compressible viscous flow theory overestimates the flow rate considerably at 
all pressure ratios except very near 1. The classical theory predicts a limiting mass 
flow of 0. 02 kg/sec (0. 040 Ibm/sec), about 250 percent higher than the measured flow 
rate. However, if the limit is chosen as the mass flow when the Mach number is 1, 
predicted flows are about 80 percent higher than observed experimentally for choked 

flow. 


Geometry Effect on Leakage 

Figure 38 shows the effect of seal-face radial length on leakage flow for film thick- 
nesses from 2. 5 to 25 pm (0. 1 to 1 mil). The four seal-face radial lengths used in the 
study were 0. 0127, 0. 051, 0. 127, and 0. 254 cm (5, 20, 50, and 100 mils). These are 
representative seal-face leakage flow lengths. From strictly a leakage point of view, 
the longest leakage path possible is most desirable. (However, the leakage path length 
must be compromised from a force-balance point of view when surface deformations oc- 
cur. ) Also shown in figure 38 is the relation of leakage flow rate to film thickness for 
a 0. 00102-cm (0. 4-mil) radial length. At film thicknesses greater than 6. 35 pm (0. 25 
mil) the leakage coincides with the leakage obtained for a knife edge by usii^ the theo- 
retical orifice flow equation. (Knife-edge flow sees no frictional resistance from the 
walls. ) Also, the linear variation of leakage flow with film thickness increases. For 
this narrow seal-face radial length, when flow behaves as knife-ec%e flow, the velocity 
entrance loss coefficient is the same as the flow discharge coefficient used in orifice 
flow analysis. 

Figure 39 shows mass leakage flow rate as a function of pressure ratio for a fixed 
gap of 10. 2 pm (0. 4 mil) and a varying seal-face radial length. The knife edge, as ex- 
pected, has the largest leakage flow rate. These values have been calculated from the 
theoretical orifice flow equation. The computer program result for a radial length of 
10. 2 pm (0. 4 mil) agrees with the theoretical orifice flow equation. As the radial length 
increases, the critical pressure ratio for choking decreases, as indicated by the sonic 
line in figure 39. Again the advantage of longer seal faces is apparent. 

Figure 40 shows the pressure profiles along the seal-face radial length found by 
usii^ the quasi-one-dimensional computer prc^ram for parallel surfaces only. For 
film thicknesses equal to or less than 5. 1 pm (0. 2 mil), the profile does not depend on 
film thickness. The flow is subsonic, as is found from the isothermal viscous flow 
analysis and shown in figure 32. When the flow is choked and the mean film thickness 
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increases, the entrance pressure loss becomes significant. (The entrance pressure 
loss is negligible in subsonic viscous flow. ) Also, under choked flow conditions, the 
exit pressure is larger than the sump pressure and increases with film thickness. As 
stated previously, the pressure decreases through expansion waves to the sump pressure 
in the outer cavity. (This does not significantly affect the axial force balance. ) 


Opening Force 

The area under the pressure gradient curves represents a force which tends to open 
the seal. This opening force is found by an integration of the pressure distribution 
across the seal-face radial length. The marked effect of seal-face distortion can be 
shown by plotting the openir^ force as a function of mean film thickness (fig. 41). For 
illustration, the opening force is plotted for five angular face distortions or (0, ±0. 001, 
and ±0. 002 rad) for the design point of 45 N/ cm^ abs (65 psia) and 211 K (100° F). As 
shown in figure 41 the effect of angular distortion is more pronounced as the film thick- 
ness decreases; the larger distortion (0. 002 rad) has a greater effect on opening force 
than the smaller distortion (0. 001 rad). As stated previously, for parallel faces the 
opening force is independent of film thickness; for convergent faces the force increases 
as the film thickness decreases (positive film stiffness); and for divergent faces the 
opening force decreases as the film thickness decreases. 

For classical viscous compressible flow, the opening force is independent of the 
fluid properties and the film thickness. This gives the well-known result that parallel 
surfaces yield no film stiffness. However, this is not the case when the classical com- 
pressible viscous flow theory is no longer valid. Figure 42 shows seal opening force as 
a fimction of film thickness obtained by using the computer prc^ram (ref. 9), for both 
isentropic and 0. 6-entrance-loss-coefficient cases. For small film thickness (corre- 
sponding to M < l/y/y), the force is constant, as predicted by classical analysis. How- 
ever, as film thickness and Mach number increase, the force actually increases slightly 
for the isentropic entrance case. This is a condition of negative film stiffness. It would 
be undesirable to operate in this region unless an auxiliary film stiffness generating 
device, such as self-acting lift pads, could make the overall seal stiffness positive. 

Also, in figure 42 the opening force attains a peak value and then sharply decreases with 
increasing film thickness for the 0. 6-entrance-loss-coefficient case. 


Center of Pressure 

Figure 43 illustrates the radial center-of-pressure variation obtained with the 
quasi-one-dimensional-flow computer program (ref. 9). At low values of film thickness 
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(subsonic flow), the center of pressure is the same as the value obtained from the iso- 
thermal viscous flow analysis (47. 37 cm; 18. 65 mils), which is valid for Mach numbers 
less than 1/^. For film thicknesses larger than the critical chdcing film thickness, 
the center of pressure increases at a lesser rate and, eventually, at larger film thick- 
nesses than shown, would be asymptotic to a dimensionless value of 0. 500. 


Continuum and Noncontinuum Flow 

Generally, gas film seals operate in the continuum flow regime, where the gas is 
very dense; that is, the molecular mean free path is very small compared to the seal 
surface separation distance (film thickness) h. The parameter relating the mean free 
path with h is called the Knudsen number. The IQiudsen number is defined as 

^ Molecular mean free path ^ Constant x T 
hj^ Mean film thickness h^j^P 

The Chapman relation expresses the mean free path as a function of the pressure 
and gas properties 


..-Ml 

2P 

The relation of molecular flow rate to volumetric flow rate can be found from 

N = ^ Q 
kT 


where 

N molecular flow rate, molecules/ sec 

k Boltzmann constant, k = 

Nq Avogadro number, 6. 027x10^^ molecules/mole 

Values of the Knudsen number indicate the regime: 

(1) Continuum flow, Kn < 0. 01 

(2) Transition (slip) flow, 0. 01 < Kn :s 1 

(3) Molecular flow, Kn > 1 

The equations just presented are for the continuum flow regime. 
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Noncontinuum flow (discussed in ref. 19) can occiu* in seals operatii^ in a vacuum 
environment. Molecular mean free paths of liquids are sufficiently short that molecular 
flow need not be considered for most practical seals. The boundary lubricating mode 
predominates over mean-free-path fluid effects. 


Flow Through Porous Media 

Porous materials are used as sealing elements. (For example, the porosity of 
carbon is highly dependent on the process, type, treatment, and impregnation. Carbon 
manufacturers usually specify the porosity. ) In high-temperature applications, porous 
materials are sometimes used for their abradability properties and/or because trans- 
piration cooling can be used. The equation that governs horizontal, isothermal flow of 
a gas is 


fx dx 

where k is the permeability in darcys (length- squared units). 

Permeability is a measure of the ease with which a fluid flows through the medium 
and is a statistical average of the fluid conductivities of the flow chaimels in the medium. 
Hence, permeability is empirically determined. 

Porosity is defined as the void volume, or the voliune of pore space divided by the 
total volume of the medium. The void volume is found by gravimetric techniques or by 
gas-e}q)ansion methods. The amount of liquid needed to saturate the dry medium is 
measured. 


INVISCID FLOW EQUATIONS 

Knowledge of orifice and nozzle flows is important in fluid sealing. Various modi- 
fied forms of these equations are used in predictii^ flow in seals, such as labyrinth 
seals, and in entrance regions and as limiting cases of maximmn theoretical leakage 
flows for a given cross-sectional area. The flow is analyzed as inviscid (neglects vis- 
cous friction forces); however, as will be seen, real flow effects such as those due to 
friction can be accounted for by various empirical coefficients. 

The inviscid flow conservation-of-momentmn equation for one-dimensional incom- 
pressible flow (eq. (136)) with the viscous friction term eliminated is 
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-A dP = M du = puA du 


(138) 


Integratii^ between any two points 1 and 2 gives 



dP = 



u du 


or 


2 2 
puj pu„ 

Pl+— Po+ — 
1 2 


This is the classical incompressible Bernouilli equation. 


(139) 


(140) 


Thin- Edge Orifice Flow 


For orifice flows, the entrance velocity is generally negligible. Thus, 


/2(Pl - 


Thus, 


^deal " -^2^2 ' H 




Po) 


(141) 


(142) 


However, in real flow situations, this theoretical value is not achieved because of 
the ’’vena contracta" effect (fig. 44(a)). This effect is accounted for by an empirically 
determined discharge coefficient Cjj which is a function of pressure ratio and Reynolds 
number. An ASME publication (ref. 20) contains data on many thin-ecfee orifice flow 
conditions. Thus, the mass flow rate is expressed as 


M = pC 


A 


- P9) 


(143) 


where 
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Actual flow rate 
Ideal flow rate 


Figure 45 shows that the orifice discharge coefficient increases with larger pres- 
sure ratios. This can be partially explained by examining figure 46 (taken from ref. 21) 
Figure 46(a) shows a rounded orifice at low Reynolds number. The air molecules are. 
able to "stream in” smoothly with a negligible vena contracta. Li contrast, the sharp 
entrance edge shown in figure 46(b) results in a stream area that is about 60 percent of 
the orifice area. However, at high Reynolds numbers (high approach velocities) the 
stream area approaches the orifice area because of the momentum of the jet air (fig. 
46(c)). 


Gravitational-Head Orifice Flow 

Orifice flow through a container with only a gravitational head driving the fluid can 
be found from Bernouilli's equation along a streamline. Refer to figure 47. 

pgh=P2=ipu2 


The exit velocity is thus 


Ug = 

This is known as Torricelli's law. The voliune flow rate is 

Qideal = 


Nozzle Flow 

Nozzle flow is generally characterized by a decrease in area between entrance and 
exit (fig. 48). The continuity equation is 

Q = A^uj = A 2 U 2 (144) 

Substituting the inlet velocity into the momentum equation, where Uj 0, yields 
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(145) 


u2-„2 = 2(Pj-P2) 

P 


In terms of the ratio of exit area to inlet area, this is 




2(Pi - P 2 ) 


(146) 


where 



Thus, 


Q = 


CpAa 2(Pi-P2> 
P 


The mass flow rate can be fovind from 


(147) 


M = pQ 

CpA2P J2(P^ - Pg) 

For compressible flow, it is common practice in flow measurements to generalize 
equation (148) by multiplying the right side by an expansion factor Y and specifying the 
specific density, p = pp as that at the inlet. Hence, both compressible and incompres- 
sible flow can be expressed by 

CpAgPiY /2(Pi - P2) 

' » Pi 

where is the incompressible flow discharge coefficient. For an incompressible 
fluid, Y = 1. It is common to theoretically determine values of Y for flow nozzles and 
venturi meters and to use experimentally determined values of Y for orifice meters. 
Values can be found in reference 20. An example of a theoretical calculation for the 
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expansion function can be found in appendix C. 

In fluid sealing, it is generally desirable to restrict or prevent leakage. Thus, a 
small discharge coefficient is desired, which is opposite to most fluid-flow applications. 
As shown in figure 45 the nozzle discharge coefficient is much higher than that of an 
orifice for the same conditions. This is due to a larger vena contracta, as illustrated 
in figure 44(b) compared with figure 44(a). Also the discharge coefficient in nozzle flow 
accounts for viscous friction of the nozzle wall as well as the vena contracta. 

Equation (139) illustrates the desirability of using a low-density gas to reduce the 
mass flow rate: 


M oc 

A gas at a low pressure (e. g. , approaching vacuum conditions) can have a low mass flow 
rate even though the gas volume flow rate may be high. This philosophy is used in some 
sealing applications where the gas is expanded to a low-pressure stage before it leaks 
out of the seal to ambient conditions. 


Flow Function Approach 


An alternate formulation of the leakage-flow equations is the flow function approach. 
This approach is extensively used in hydraulics and in sealing systems, especially for 
lab 3 rrinth seal flows. This formulation is based on the one-dimensional inviscid flow 
assumption but accounts for viscous and other effects by coefficients which are usually 
empirically determined. 

For an orifice flow, equation (149) becomes (setting /3 = 0) 


2(Pj - P ) 
M = Cj^Ap^J i 2_ 

* ^1 


For a perfect gas, equation (150) becomes 


LM Piy 


(150) 


(151) 


or 
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(152) 



We define the flow function as 


APfl 


The flow function is one dimensional and is a function of the pressure ratio and the par- 
ticular overall geometric configuration. The pressurized fluid conditions are now des- 
ignated by the subscript 0. The flow function is a measure of the flow effectiveness. 
The dimensionless mass flow equation (152) can be rewritten in the following form: 



= aycpC 


where 


G = M/A a convenient parameter for sizing orifices and knife edges 


a - Cd 
r 

c 


discharge coefficient 


recovery factor, used for straight-through staged restrictions (e. g. , in 
stages of orifices or labyrinth knife edges) 


conversion constant 


(e.g.,yS) 


The famous Egli formula (ref. 22) used in estimating labyrinth seal leakage is of such 
a form. The functions a, y, and ip from that formula are plotted in figure 49. 

The industrial gas compressor industry uses a formula 


M = KCj^ApUj^ 

where the coefficient K is empirically determined. 


Labyrinth Seals 


High friction losses due to turbulence are produced in throttling or restrictor seals 
such as labyrinths. The sealed gas expands as it passes through orifices or knife edges 
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formed by successive restrictions. An energy conversion process takes place. The 
stored gas (internal) energy (due to pressure and temperature) is converted into kinetic 
energy in the flow process. (The energy contained in the gas by its temperature and 
pressure state is called the gas enthalpy. ) By designing the geometry within a given 
configuration, the kinetic energy is converted into heat; otherwise, a near full recovery 
would occur. A common way of dissipating the kinetic energy is by changing the flow 
direction between successive restrictions. A ’’stepped" labyrinth illustrates the use of 
this idea. 


ENTRANCE FLOWS AND LOSSES 

It is common in lubrication theory to neglect the entrance region entirely. This is 
reasonable for the slow viscous flows which characterize lubrication flows. Brodoia and 
Osterle (ref. 23) performed a numerical analysis of plane Poiseuille developing flows. 
They found that the entrance lei^th varied approximately as (fig. 50) 

Xj^^O. 044 hRe (153) 

Thus, it is seen that for low Reynolds number flows this region is negligible. Kawaguchi 
(ref. 24) studied entrance losses of flow between parallel disks and concluded that they 
were the same for low Reynolds numbers as for laminar flow. However, there are con- 
ditions when entrance losses are substantial, namely, when the flow is high Reynolds 
number (turbulent) and when the flow becomes choked in compressible flow. 

There are several causes for pressure drops in the face-seal entrance region. An 
abrupt geometric change at the seal passage entrance can result in a nonuniform profile 
because of flow turning, separated flow, and the vena contracta effect. This nonuniform 
profile will cause an entrance pressure drop (fig, 51). Kawaguchi (ref. 24) found that 
the entrance loss for turbulent flow depends strongly on the ratio of film thickness to 
inlet radius and on the radius of inlet curvature but depends very little on the entrance 
Reynolds number. The pressure drop in the flow-development length is higher than that 
in the fully developed flow region because of two effects (ref. 25). The first effect is 
higher wall shear caused by higher transverse velocity gradients. The second effect is 
the momentum increase as the velocity distribution becomes less uniform. These en- 
trance pressure losses can be accounted for by using entrance loss coefficients. 

The role of the entrance velocity loss coefficient may be better understood by con- 
sidering the incompressible Bernouilli equation relating the stagnation reservoir pres- 
sure with the static and dynamic pressure at the seal entrance: 
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(154) 


^0 " ideal 


The actual velocity at the seal entrance is less than ideal; it can be expressed as 


"l - ^l'^I, ideal 


(155) 


where is the entrance velocity loss coefficient (ref. 26) and is experimentally de- 
termined. Thus, the entrance pressure can be found from 



(156) 


The head loss coefficient k’ is commonly used in hydraulics and reported in the 
literature The velocity loss coefficient can readily be related to the head loss coeffi- 
cient by the following equation: 


C, =-4= 


(The head loss coefficient is sometimes also called the flow resistance coefficient. ) 
Fleming and Sparrow have shown in reference 25 that the bulk of the entrance losses for 
duct flow occur in a small region very close to the duct inlet. Laminar incompressible 
entrance loss coefficients have been calculated and measured for a variety of duct shapes 
(e. g. , refs. 25 and 27). Turbulent loss coefficients are not so widely reported; how- 
ever, they may be calculated for parallel-plate channels from the information in refer- 
ence 28. The values for head loss coefficients k' for turbulent flow are generally less 
than 20 percent of those for laminar flow. See the qualitative comparison shown in fig- 
ure 50. Incompressible loss coefficients may be adequate for use with compressible 
flow (ref. 11), 

The head loss approach is also used to analyze sealing systems. The head loss is 

H = ^ 

P 
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Thus, the inviscid flow velocity equation becomes 


u = 


where Cj^ is the entrance loss coefficient and 




The entrance loss coefficient is related to the head loss coefficient (or flow resistance 
coefficient) by 


■ 


Vk’ + 1 


Hj =— (k’ + l) 


For viscous effects, equation (88) 


dx 


can be integrated to the form 


' ^2 u^ 4fL ^ jj 

p 2D 


Thus, the total head loss for both entrance and viscous losses is 


H = /k’ + 1 + 


where 


J 
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Here 


K = 


1 




(157a) 


A head loss equation for labyrinth seals that is especially useful for seals with 
honeycomb-cell rub strips is 


H = 



where n is the number of throttling stages. 


FLOW REGIMES 

As mentioned in the previous sections, the appropriate equations for seal analysis 
depend on the flow regime. There are many possible flow regimes because of the wide 
spectrum of seal operation. The most common ones will be discussed in this section. 

The leakage flow regimes are defined by either the appropriate Reynolds number or, 
in the case of a gas, the molecular mean free path ler^h. The flow regime must be 
established in order to determine the applicable theoretical leakage relation. Since the 
leakage flow is required for the Reynolds number calculation, it is necessary to assume 
a flow regime and to iterate the leakage and Reynolds number calculations. 

Consider flow in a seal passage caused by a pressure differential. When the flow 
is laminar, the fluid flow is characterized by laminas, or layers (fig. 52). This flow is 
predictable analytically, as we have seen for flow in long passages. The chief charac- 
teristics of this viscous laminar flow are that the velocity is independent of the internal 
surface roughness and that the velocity distribution across the gap is parabolic in form 
with the maximum velocity at the mi^ap position. At flow velocities above a critical 
condition, the flow pattern changes to one of turbulence characterized by eddies and vor- 
tices (fig. 53). The energy losses now become affected by the surface roughness, and 

flow prediction requires the use of experimental data. Turbulent flow is characterized 
by velocity components which have random turbulent fluctuations imposed upon their 
mean values. A stream of dye or ink inserted in a laminar flow will streak out a thin 
line and will always be composed of the same fluid particles. However, in turbulent flow 
the dye line would quickly become tangled and mixed in with the fluid as it flows along, 
and we would see myriads of threads and clouds widening and dispersing as the fluid 
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flows along. Figure 54 qualitatively compares laminar and mean turbulent velocity dis- 
tributions. The turbulent mean velocity profile is more elliptic (blunt nose) and is char- 
acterized by a much higher shear stress (due to the slope of the velocity profile at the 
wall). Further details can be found in any fundamental fluid mechanics textbook. 

Figure 55 qualitatively shows how the pressure drop per unit flow length varies as 
the fluid velocity increases. There is a sharp increase in pressure drop when turbulent 
flow occurs. This sharp increase in pressure drop, which is a measure of the shear 
force, is characteristic of turbulent flow. 

The exact transition point for any particular set of conditions is uncertain, but the 
limits of the laminar-turbulent flow transition regime can be defined by the value of the 
nondimensional parameter known as Reynolds munber Re, where 

Re=£HD 

M 

The minimiun transition Reynolds number for pressure flow is 2300. This has been ob- 
served to be the case for many geometries and experimental conditions in fluid mechan- 
ics and hydraulics. For flow between narrowly spaced plates and coaxial annular disks 
the hydraulic diameter is D = 2h. 

For highly viscous liquids, such as oil, used in most sealing applications the flow 
will be laminar. Lower viscosity liquids, such as water and especially liquid metals, 
can operate in the turbulent flow regime at very high speeds. For gases, such compres- 
sibility effects as fluid choking appear to be important before the flow becomes turbu- 
lent. Hence, the Mach number, as well as the Reynolds number, will be important. 

Shear flow . - H there is relative rotation of the seal surfaces, as in mechanical 
face seals, turbulence can be induced by this rotation. The rotational flow component of 
the velocity may be responsible for the flow becoming turbulent in the circximferential 
direction (shear flow direction), which means the entire flov^field (radial pressure flow) 
will be turbulent. 

For narrow gaps (close clearances, where the boundary layers are merged), the 
critical rotational Reynolds number for transition appears to be the simple Couette flow 
transition Reynolds number. Thus, Re^ = rfih/i/=^ 1900 for transition (ref. 29); the 
smaller the gap, the longer the transition speed is delayed. Therefore, the wall has a 
stabilizing effect on the flow, and this should be considered in selecting the design gap 
(e. g. , for a 30-cm- (1-ft-) diam air seal operating at room temperature, the transition 
speed would be 71 500 rpm for a 25-jim (1-mil) gap but would be 179 000 rpm for a 
10-pm (0. 4-mil) gap). 

Some seals, such as vaneless slinger seals, operate at larger clearances in order 
to seal high pressures at high speeds. (Heat dissipation must be minimized, so clear- 
ances are relatively larger than with face seals. ) 
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Daily and Nece (ref. 30) performed experiments on enclosed rotating disks. They 
identified four basic flow regimes: 

(1) Regime I; laminar flow, close clearance 

(2) Regime II: laminar flow, separate boundary layers 

(3) Regime III: turbulent flow, close clearance 

(4) Regime IV: turbulent flow, separate boundary layers 

Similar experiments were performed by Ketola and McGrew (ref. 31), who identified 
the flow regimes for a partially ei^aged seal. A schematic of their seal is shown in fig- 
ure 56. The approximate flow regimes are shown in figure 57. The characteristic 
lei^h in the Reynolds number is the radius r instead of the gap h. This is the char- 
acteristic length for a rotating disk in an infinite fluid medium. For a highly polished 
disk the transition Reynolds number 


2 

3x10^ 

V 


Since transition from laminar to turbulent flow depends on the disk radius, part of the 
seal (disk center region) will always be operatii^ in the laminar flow regime. Whether 
or not the rotating boundary layers are merged can be estimated by applying von Kar- 
ma's formula (ref. 32). For laminar flow, 

6^2. 58 

1S2 

Theoretical radial and circumferential velocity profiles are shown in figure 58 for 
various Reynolds numbers. These are taken from Lance and Rogers (ref. 33). For the 
boundary layer case (Re = 441), theoretical solutions are compared with experimental 

results in figure 59. 

Annular flow between cylinders . - When inner and outer cylinders are very closely 
spaced, the flow remains laminar until the rotational transition Reynolds number is ex- 
ceeded. However, if the gap is sufficiently large, a secondary flow regime called Tay- 
lor vortex flow can occur when the critical speed of rotation is reached. These vortices 
are ring^aped around a rotationary axis, about twice the gap size thick. Taylor vor- 
tex flow is illustrated in figure 60, and the vortices (paint pigment particles immersed 
in oil) are shown in figure 61 (ref. 34). Taylor (ref. 35) analyzed mathematically the 
critical speed of rotation at which these vortices are formed, and he confirmed it by his 
experiment. The criterion is that Taylor vortices occur when the Taylor number 
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Ta = 







Rj + R2y 


> 41.2 


Kaye and Elgar (ref. 36) observed that Taylor vortices keep their shape until the 
speed of rotation is increased from a critical value to a certain value. Taylor vortex 
formation has been observed to be entirely different from turbulent flow formation. 

When the speed of rotation of the inner cylinder is increased by a large amount, the flow 
in the annulus becomes turbulent and vortices are superimposed on the turbulent fluc- 
tuations. When the speed of rotation is increased further, the Taylor vortices are ob- 
served to first distort and then completely disappear, and a purely turbulent flow re- 
sults. 

Kaye and Elgar (ref. 36) identified four modes of flow between coaxial rotating cyl- 
inders when there is a net axial flow: 

(1) Purely laminar 

(2) Purely turbulent 

(3) Laminar plus Taylor vortex 

(4) Turbulent plus Taylor vortex 

These flow modes depend on the inner cylinder rotational speed and mean axial flow ve- 
locity. A schematic representation is shown in figure 62. 

Pressure flow with rotational flow . - When there is both circumferential shear flow 
and radial pressure flow present (as in a face seal), the net flow is a spiral flow. In 
the laminar flow regime the shear flow and the pressure flow can generally be treated 
independently. This was concluded from experiments performed by Yamada (ref. 37), 
who studied axial flow through an annulus with the inner cylinder rotatii^. The rota- 
tional flow is always important for determining the power loss due to viscous shearing 
and for determining the transition to turbulent flow. Of course, it can be responsible 
for rotationally induced turbulence. As previously seen, the centrifugal component in 
the radial direction can contribute to the seal leakage, especially for liquids. Yamada 
also discusses various effects of rotation and the effect of the axial flow becoming turbu- 
lent. 

Figure 63 illustrates a simplified envelope of possible seal-face compressible flow 
regimes when there is circumferential rotation as well as radial leakage flow (ref. 38). 
As previously stated, the Mach number is the radial flow parameter, and a simple 
Couette flow transition Reynolds number is used as the rotational flow parameter for 
very narrowly spaced disks representative of face seals. Theodorsen and Reigier 
(ref. 32) found that laminar-turbulent transition and drag were independent of compres- 
sibility for a disk rotating in an infinite medium. Although not verified, it is speculated 
that this is also the case in face seals. 
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MECHANISM OF FILM PRESSURE GENERATION BETWEEN 
RELATIVELY MOVING SURFACES 

Full-fluid-film lubrication implies the existence of an interposed full film of fluid 
that physically separates the mating surfaces of the two members which traditionally 
comprise a bearing. However, fluid-film lubrication plays a key role in sealing as well. 

Three generally recognized modes of lubrication are categorized as the sliding- 
surface bearing type (ref. 39): 

(1) Hydrodynamic - in which fluid-film, load-carrying capacity is achieved by mo- 

tion of the bearing members in the plane of the mating bearing surfaces 

(2) Hydrostatic - in which fluid-film, load- carrying capacity is achieved by the in- 

troduction of externally pressurized flow 

(3) Squeeze film - in which fluid-film, load-carrying capacity is achieved by motion 

of the bearing members normal to the plane of the mating bearing surfaces 
These three modes of fluid-film lubrication are illustrated diagrammatically in figures 
64(a) to (c). A possible fourth mode, termed "hybrid operation, " is illustrated in fig- 
ure 64(d) and derives its fluid-film, load-carrying capacity by a combination of two or 
more of the previously mentioned modes of lubrication. 


Hydrodynamic Pressure Development 


The classical lubrication equation is called the Reynolds equation. It is the differ- 
ential equation which relates viscosity, film thickness, pressure, and sliding velocity. 
This equation is derived in many lubrication textbooks. In the section Flow Between 
Parallel Plates in Relative Motion, we have already used the same assumptions that are 
used in deriving the Reynolds equation. Appendix D contains a brief description of 
Reynolds equation taken from Bisson and Anderson (ref. 40). 

The Reynolds equation is essentially a flow balance equation. The viscous flow 
assumption (that the viscous friction force is balanced with the pressure gradient force) 
is made. The momentum equation is then substituted into the mass conservation equa- 
tion and integrated across the gap. The two-dimensional Reynolds lubrication equation 
is 


X- Direction z- Direction 

pressure pressure 

flow flow 
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to wall 
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Physical explanation . - A physical understanding of pressure development in an oil 
film can be achieved by examining the pressure conditions that must prevail to maintain 
flow continuity with various film configurations. 

Consider two parallel plane surfaces separated by a lubricant film. As shown in 
figime 65, the lower surface is stationary and the upper surface is moving with a ve- 
locity U in a direction parallel with the surfaces. The local flow velocity in the film 
u varies uniformly from zero at the stationary surface to U at the moving surface. If 
the plates are very wide in the direction perpendicular to the direction of motion, the 
rate of flow that crosses the left boundary AB will be equal to the flow that crosses the 
right boundary CD. The flows that cross these two boundaries result only from velocity 
gradients. Since they are equal, the continuity requirements are satisfied without any 
pressure buildup. Therefore, no pressure buildup wiU result when two parallel surfaces 
are moving relative to each other only in a direction parallel to the surface. This was 
shown in the section Flow Between Parallel Plates in Relative Motion. 

Next, consider the two plane surfaces inclined slightly with respect to each other, 
with one of the surfaces in motion in a direction parallel to its surface (fig. 66(a)). 

Again, as with parallel surfaces, the local flow velocity u in the film varies uniformly 
from zero at the stationary surface to U at the moving surface. Assiune that the plates 
are wide in the direction perpendicular to the direction of motion so that no flow occurs 
except in the direction of motion. Then the velocity flow that crosses the boundary AB 
will be greater than the flow that crosses the boundary CD. The requirements of flow 
continuity thus are not satisfied since no sources or sinks exist between the two bound- 
aries. It can be surmised that there will be a pressure buildup in the oil film until the 
continuity of flow is satisfied. 

The flow which results from the pressure buildup is sketched schematically in fig- 
ure 66(b). The local flow velocity that results from pressure gradients is zero at both 
surfaces since no interfacial slip is assumed. The velocity is a maximum at a point 
halfway between the two surfaces. The velocity profiles are parabolic. Because pres- 
sure flow is always from a region of higher pressure to a region of lower pressure, the 
flow is out of the bearing at both boundaries. 

The pressure flow at the boundary AB opposes the velocity flow, while the pressure 
flow at CD is in the same direction as the velocity flow. If the velocity and pressure 
flows are superimposed, the flow profiles that result are as shown in figure 66(c). The 
required equality of flow across AB and CD can now be maintained with the proper com- 
bination of velocity and pressure within the film. At some section EF within the film 
there is no pressure flow, and the velocity profile across the film is again linear. This 
is the point of maximum pressure or the point at which there is no pressure gradient. 

It was deduced previously that some inclination of two plane surfaces in relative 
motion is required in order to achieve a pressure buildup. This deduction can be qual- 
ified further by requiring that the film be convergent in the direction of motion. A di- 
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vergent film will produce not a pressure buildup (eq. (60)) but a subambient pressure 
reduction. A convergent film in a hydrodynamic bearing is called a wedge film. A 
wedge film is one of two ways in which pressure or, equivalently, load capacity can be 
generated in a hydrodynamic bearing. 

Narrow -seal approximation. - The narrow-seal approximation is commonly used in 
analyzing hydrodynamic face seals. This approximation neglects the circumferential 
pressure gradient contributions to the circumferential mass flow rate. This is analo- 
gous to Ocvirk and DuBois' short-journal-bearing theory (ref. 40): 
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This assumption amounts to claiming that the majority of the circumferential mass flow 
is caused by tangential viscous shear (accounted for by the right side of the equation). 

For narrow seals with constant boundary pressure around the perimeter, it is usually 
unlikely that large circumferential pressure gradients can exist within the seal except 
possibly in narrow interior regions along the boundaries. 

Step bearing. - A common configuration that generates hydrodynamic pressure is the 
step bearing configuration shown in figure 67. The main assumptions are that hj « 
and h 2 « L 2 . The flow is assumed to be viscous, two-dimensional laminar flow. At- 
tach a coordinate system to the slider, as shown in figure 67. Let the stationary bearii^ 
surface move with a velocity U in the negative x-direction. Hence, the boundary con- 
ditions are 


u = 0 at z = 0 

u=-U atz = h 


and the gage pressure is 


P= 0 

at the inlet X 2 = L 2 and the exit Xj^ = 0. The pressure at the step Xj = or X 2 = 0 
is unkn own. This pressure must be found by invoking the continuity condition that 

Ql = Q2* 

Assuming that general Couette flow governs the flow along flow passages Lj^ and 
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L 2 and using equation (87) give the continuity condition as 





12ju L. 


Uhj 

2 



12/xLp 


2 


The pressure in the land and the groove vary linearly with x, and a triangular pressure 
results. The peak pressure at the step Pg is found by equating Qj and Q 2 . The re- 
sult is 


P 


s 


6(uU(h2 - hj) 



When this pressure is known, the entire distribution across the land and groove are 
known. 

The total hydrodynamic load generated can be found by integrating the pressure dis- 
tribution curve. The result is 


F = — § W(L^ + L2) 

2 

If the groove depth and film thickness were changed, the slider would rise or sink as the 
load was decreased or increased, respectively. Hydrostatic step seals and viscoseals 
use this life mechanism in their operation. 

Compressibility effects . - Hydrodynamic lift concepts apply to gases as well as to 
liquids; however, compressibility effects (ref. 41) must be accounted for. The same 
assumptions apply as were stated for classical viscous compressible flow. Namely, 
isothermal, laminar, viscous flow of a continuous fluid with negligible inertia is as- 
sumed. 

Since gas lubricant films are very thin and in contact with relatively large masses 
of metal, the films can be considered to be isothermal. Therefore, changes in local 
density in the film result primarily from changes in the local pressure as predicted by 
the perfect -gas law. This Interdependence alters the pressure distribution from the 
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shape it would have if it were incompressible (fig. 68(a)), Consequently, the degree of 
gas compression affects the bearing performance. 

The compressibility number A is a measure of the density variation importance. 

In the presentation for compressible flow the Mach number was stated to be the com- 
pressibility effect criterion. This parameter was an indication of when the inertia ef- 
fects were important. In the viscous laminar flow regime, inertia effects are negligible 
and the Mach number is not used. Also the compressibility parameter can be inter- 
preted as 

A-^ 

P 

For small A, the hydrodynamic effect is the same as for incompressible fluids. 
When A is large, figures 68(b) and (c) show that incompressible theory does not apply 
and that compressibility causes a deterioration in performance. Figure 68(c) illustrates 
a barrier which is imposed by compressibility on most types of gas bearings. The load 
capacity of the gas bearing does not indefinitely increase with increasing speed but 
levels off. In contrast, the load capacity of the liquid bearing is directly proportional 
to the speed. 

Although this brief discussion applies to a basic fixed pad, the same effects are 
present in most other hydrodynamic bearing configurations. This includes spiral 
groove, herringbone, and tilted pad bearings. Further discussion is beyond the scope 
of this presentation. The solution procedures are similar to those already mentioned; 
however, compressibility does complicate the solution. For more information, refer 
to the textbook by Gross (ref. 15) or to reference 41, 

Since gas viscosities are of the order of 1/ 1000th that of oil, hydrodynamic lift is 
much less for gases. The lift is also reduced because of compressibility effects. That 
is, as the pressure rises through the bearing, the density increases. The changing 
volume flow through the bearing modifies the pressure characteristic, giving lower load- 
carr 3 ^ng capacity. If the compression ratio Pjj^ax^^a close to 1, compressibility 
effects become insignificant and performance efficiency is imposed. 

Squeeze-Film Concepts 

Relative seal surface vibrations are inherent in all rotating seal systems. A typ- 
ical vibration is one surface oscillated at a high frequency perpendicular to the plane of 
the surfaces. This phenomenon can best be investigated by studying a liquid in a narrow 
gap where one surface is oscillated at a high frequency. Figure 69 illustrates two par- 
allel plane surfaces which are separated by a lubricant film, with the upper surface 
moving toward the lower. Lubricant is being squeezed out of the volume between the 
surfaces. Since the surfaces have no velocity in the direction of flow, the only cause of 
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this flow is a pressure gradient. Relative motion of two surfaces separated by a lubri- 
cant film toward each other will then result in a pressure buildup, as shown in figure 69. 
In a hydrodynamic bearing this film is called the squeeze film. 

Hydrostatic cushioning . - Squeeze-film concepts play an important role in impact re- 
sistance. Generally, good impact loading resistance is obtained if the duration of im- 
pact is less than the squeeze-film duration. A viscous lubricant cannot be instantane- 
ously squeezed from between two surfaces which are approaching each other. It takes 
time for these surfaces to meet. During that interval, because of the resistance of the 
fluid (lubricant) to extrusion, a pressure is built up and the load is actually supported by 
the fluid film. This is called hydrostatic cushioning, a special case of the squeeze-film 
concept. If the load is applied for a short enough period, the two surfaces may not con- 
tact. The governing form of the Reynolds equation for this case is 

li. /(.SrpiP 

r 9r y 9r 


12m 
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(158) 


is 


For one circular disk approaching another, the flow rate for an incompressible fluid 

q = TTT^V 


where v is the relative velocity of the disk approach. K we assume laminar, viscous 
flow of the fluid effluxing out the boundary. 
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Applying the boundary condition that at r = R, P = 0 yields 

P = 3^(r2 
u3 
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This equation is for a paraboloid, for which 


Pav = 


max 


and P = when r = 0. Thus, 


P 
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_ 3 uvR^ 
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The squeeze-film load can be found from 


W - ttR^P = 

av o 

2h'^ 


Observe the functional relations of the variables. Notice the inverse cubic dependence 
on film thickness. 

The squeeze“film time can be found from the load equation after substituting 

dh 

V = 

dt 


Hence, 
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An example case (castor oil) is 


hj = 0. 0004064 cm (0. 00016 in. ) p = 2. 43xl0‘3 N-sec/in. ^ (5. 08x10'^ Ibf-sec/in. 

h 2 = 0. 0001270 cm (0. 00005 in. ) W = 0. 885 N (0. 199 Ibf) 

The calculated s<jueeze“film time is 71. 6 minutes. This has been verified experimen- 
tally. Further discussion of hydrostatic cushioning can be found in references 41 and 42. 

Relative sinusoidal vibration of surfaces . - Since the disk vibrates sinusoidally at a 
frequency cu. 


h = + e cos wt ( 159 ) 

where e is one-half the amplitude of the movement about a mean film thickness h 
Equation (159) can be written as 

h = hQ(l + e cos wt) (160) 

where e is called the excursion ratio, 


e 


e 


h 


0 


To nondimensionalize, let 


T - Ojt 



Thus, equation (158) becomes 
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h!! _L 4*r* ^]=a 

r* dr* \ 9r*/ 


(161) 


a(p*h*) 

9t 


where (t = 12R^/iw/h^P^ is the squeeze number. 

Salbu (ref. 43) has solved this equation for the pressure distribution and the gener- 
squeeze”film forces. Figure 70 shows the dimensionless squeeze“film force for 
two excursion ratios (e = 0. 1 and e = 0. 5) and a squeeze number of 10. For the small 
excursion ratio (e = 0. 1), the variation of the squeeze-film force with time is almost 
symmetrical, which results in a zero net load when integrated over a cycle of oscilla- 
tion. At the larger excursion ratio (e - 0. 5), the variation is asymmetrical, which gives 
a corresponding load. Thus, the behavior can be summarized in the following way: K 
the excursion ratio is small, the fluid is compressed very little and the pressures are 
almost symmetrical; hence, the net load is zero. However, when the excursion ratio 
is large, on the downstroke the air is compressed and resists the movement strongly. 

On the return and upstroke the gas becomes rarefied and offers less resistance. The 
squeeze-film force is thus highly asymmetrical, and a net load capacity exists over a 

cycle. 

This behavior can be thought of in another way. At a low excursion ratio and/ or 
low frequency the viscous forces predominate over the compressibility forces. Hence, 
gas moves radially to compensate for the normal motion. On the other hand, for a high 
excursion ratio and/or high frequency the gas moves very little radially, being com- 
pressed and rarefied. Under this condition the movement is very similar to a piston in 
a closed cylinder. Hence, Boyle’s law, PV = Constant, applies. 

Consider a simple, squeeze-film, thrust bearii^. For high frequency, the film 
compresses and expands so quickly that there is very little leakage. If the very thin 
boimdary layers near the boundaries are neglected, the P and V satisfy Boyle's law 
(PV = Constant). Hence, 


P^(V + AV) = PqV = P_(V - AV) = Constant 


or 


where 


P = 



= P 1 



P 


+ 


1 


+ 


AV 

V 
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and 



V 


Hence, the average pressure due to expansion and contraction 2 AV is 



1 - 



> Pn 


and 



Thus, 


^av " ^0 



This indicates a net pressure buildup to oscillation, which causes change in the volume 
AV. The buildup in mean pressure is a second-order effect in AV, as verified from the 
more sophisticated theory. 

The net squeeze-film load can be estimated from (ref. 44) 
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w = 



where 1.15. 

A net squeeze-film force exists for an isothermal gas. However, for a liquid (in- 
compressible fluid), this is theoretically not the case. Since the sonic velocity is infi- 
nite, any disturbance of the fluid is instantaneously felt everywhere. Hence, theoret- 
ically the liquid will always move radially to compensate for the normal motion. 

However, squeeze-film bearings have operated successfully with alcohol and water. The 
reason for a net load is attributed to degasification and/ or cavitation of the fluid - a 
characteristic of real liquids. With a two-phase system, the lowest pressure is not 
reached, because of the phase change or greater concentration of dissolved gas (fig. 71). 
Hence, a net load appears. This effect probably occurs in seals and may be a mechan- 
ism for obtaining seal-face separation. 

Cavitation ajid Degasification 

As seen in the section Flow Between Parallel Plates in Relative Motion, for viscous 
flow between apparently parallel surfaces, theory predicts no hydrodynamic film. This 
would indicate rubbing contact. With the formation of a film a load-carrying capability 
is developed, resulting in a positive separation of the surfaces. Lubricant cavitation 
has been offered as a possible cause for the generation of this hydrodynamic force. 

The common explanation of the cavitation phenomenon is that it occurs when the dy- 
namic conditions at some point in a flowing fluid produce a pressure lower than its vapor 
pressure. For some liquids, however, the vapor pressures are so low that this explana- 
tion of cavitation must be modified. "Cavitation" then results from dissolved gases 
separating from a solution, which is sometimes called "degasification. " Experiments 
show that the amount of gas dissolved increases with increasing pressure and that the 
gas is released from the liquid when the pressure decreases. These situations are es- 
pecially significant for seals as most oils have low vapor pressures. Hence, degasifi- 
cation will be seen to play a key role in generating a hydrodynamic lift force; however, 
the appearance of bubbles alters the performance of the lubricant as a whole by trans- 
forming it to a two-phase fluid. This, in turn, lowers the ability of the fluid to withstand 
shear stress and to transfer heat effectively. Degradation of the lubricant results. 

Hence, there is a trade-off. The exact role of cavitation and degasification is still not 
well understood. 
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Pressure Development due to Microroughness 


Surface roughness on the microscale may be an important mechanism for generating 
an opening force when the film thickness is of the order of the surface irregularities. 

The surface irregularities are assumed to cause a pressure gradient in the oil film in 
the same manner as a pressure gradient is established across a tapered-shoe thrust 
bearing. The interface of the seal is assumed to have the form shown in figure 72. This 
is a simplification of the actual roughness pattern. Application of a simplified form of 
the Re 3 molds equation 


A 43 

dx \ dx/ 


+ 6pU— = 0 
dx 


gives the pressure distribution across one of these microsurface wedges. 

The sealing surfaces can be idealized into a smooth movii^ surface and a stationary 
surface with a roughness pattern, as shown in figure 73. Davis (ref. 45) states that if 
an idealized seal surface consists of a number of units B’AB (fig. 74) which are rigid, 
the pressure distribution would be skewed symmetrically about A and, if cavitation ef- 
fects are neglected, the net lift would be zero. However he states that if the surface is 
elastic and deforms with pressure, a new form results (fig. 74) and the pressure dis- 
tribution becomes asymmetrical. The surface now experiences a positive lift expressed 
by 


w = 


pUrf2ffRL 



where L is the seal radial length, h^ is the mean clearance, R is mean seal radius, 
r is the maximum surface separation, and f is the convergent surface spacing (fig. 73). 
Davis makes an empirical allowance for the diverging clearance by supposing the lift is 
generated by the converging gap only. 

Anno, Walowit, and Allen (ref. 46), in their experimental investigation of face seals 
with photoengraved microasperities on one of the faces, observed that a hydrodynamic 
film with small cavities was established between the seal surfaces (fig. 75). The small 
cavities began at the downstream side of each asperity and ended at the adjacent up- 
stream asperity. As the formation of the cavities prevents the generation of a hydrody- 
namic pressure less than cavity pressure, the effect of the cavities on the face- seal flow 
lies in truncatii^ the negative pressure buildup that otherwise would have balanced the 
positive pressure buildup in the vicinity of the asperity if the cavities were absent. This 
gives rise to an asymmetrical pressure distribution whose integrated effect is the de- 
velopment of a hydrod 5 mamic force. Fluids that have low vapor pressures at atmos- 
pheric conditions will generally degas the entrained vapor; water may also cavitate. 
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Pressure Generation due to ’•Wavy” Surfaces 

Hydrodynamic face separation can also be achieved from ’’wavy” seal surface ef- 
fects in radial face seals. From lapping of carbon, typically, two or three peaks and 
valleys are present. The short -bearing form of the Reynolds equation is valid and in 
polar coordinate form is 



By assuming that cavitation occurs where the pressure is subambient and the surfaces 
are diverging, a net hydrodynamic lift force is calculated. This process has been 
studied and observed by Stanghan-Batch and Iny (ref. 47) and is illustrated in figure 76. 
Similar pressure generation can be obtained by diametrically misalined seal faces. 

Other methods of pressure generation appear in the literature on seals. 

BOUNDARY LUBRICATING REGIMES 

Many seals rub during operation or during startup and shutdown. The boundary 
lubricating regime entails rubbing of seal surfaces. Before this regime is discussed, 
it is advantageous to study some fvmdamental friction and wear concepts. 

Fundamental Friction Concepts 

Several basic friction and wear concepts that have been accepted by many investi- 
gators are important in understanding the phenomena that occur in lubrication of surfaces 
in solid contact (boundary lubrication). All surfaces have irregularities (or roughness), 
however, our common measures of roughness may only define whether a surface is 
rough on a macroscopic, a microscopic, or an atomic scale. These concepts are based 
on a microscopic consideration of surface character and usually are essentially inde- 
pendent of the gross geometry of the mechanical parts. References 48 to 50 discuss 
these concepts in detail. 

A primary concept is that surface contact occurs between asperities (or protuber- 
ances). It follows that the real area of contact is a minute fraction of the apparent area 
of contact. (The apparent area of contact is defined as the gross geometry of the mating 
parts. ) As surfaces come into solid contact under light load, very few of the asperities, 
or high points, may actually touch. With increasing loads, the initial points of contact 
deform elastically and then plastically, and additional asperities assume a portion of the 
total load. The local stress of the initial asperities in contact is then at the elastic limit 
(or yield stress in compression) of the weaker of the materials in contact. Thus, the 
maximum real surface stress will usually be the yield stress, and increased total load 
simply increases the real area of contact. 
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The common laws of friction (Amonton's laws) are usually valid because plastic 
formation determines the real area of contact. However, they are not valid when defor- 
mation is primarily elastic or when the load, materials, and geometry are such that the 
real and apparent areas of contact are equal. The latter case gives rise to gross plas- 
tic flow and results in a "plowing" friction force. 

The most widely accepted expression for these concepts has been advanced by Bow- 
den and Tabor (ref. 50). This expression shows that friction coefficient tj, which is 
defined as the ratio of friction force to normal load W, is equal tojhe shear 
strength S of the softer metal in contact divided by the yield strength P of the weaker 
metal in compression. Figure 77 illustrates solid surfaces in contact. Where is 
the real area of contact, these simple concepts may be related as follows: 


= = V 


w w 




P 


Film transfer and reactive coatings . - Thin surface films, such as lubricants or 
other deposited or reactive coatings, reduce friction by lowerii^ effective shear 
strength. The important bulk yield strength factor, which determines shear area, is a 
function only of the base metal substrate so long as the film is very thin. Among the 
most common surface films important to lubrication are the surface oxides. The type, 
thickness, and character of such films are a function of the thermal history, chemical 
kinetics, and thermodynamics of the system; surface temperatxire is especially impor- 
tant. Figure 78 illustrates the influence of contamination on the friction coefficients for 
various material combinations. 

Nonmetals, such as carbon, can transfer to the matii^ metal surface. This occurs 
when a residual surface oxide is present on the metal surface. The carbon bonds to the 
oxygen of the metal oxide, and transfer occurs during sliding. The end result is carbon 
s liHing on a carbon film and a very desirable compatible situation of low friction and 
wear exists. (A mechanical lead pencil is a simple example of this transfer film. ) 

The "wearing in" of a carbon face seal can be thought of as the process in which the 
carbon transfers a thin film to the mating hard surface. Also a smearing or smoothing 
of the carbon surface microroughness may be occurrii^. However, if the mating sur- 
face is composed of a metal oxide which is thermodynamically imstable and dissociates, 
the metal can transfer to the carbon surface. This is generally a function of the temper- 
ature and load. The result is metal sliding on a metal film, an undesirable situation. 
This may be the mechanism involved in cases where carbon does not wear but the mating 
harder metal surface does. 

The energy consumed in friction during sliding of solid surfaces must be dissipated 
as frictional heat. When the real area of contact is a small fraction of the total surface 
area, the heat flux at the surface asperities may be extremely high. It has been demon- 
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strated that surface temperatures of 500° to 1000° C (932° to 1832° F) above the ambient 
temperature can easily be experienced. With chemical activity a function of tempera- 
ture, it is apparent that surface reactions will occur preferentially at the local high- 
temperature areas on the surface asperities. Excessive reactions result in chemical 
attrition or corrosive wear, so that an optimum chemical activity of the surface must be 
sought to obtain effective surface protection from both adhesive and corrosive wear. 

Reactive coatings play an important role. Many sealed fluids act as lubricants. A 
typical solid surface in contact with a lubricant is shown in figure 79. Most solid sur- 
faces in an air atmosphere contain films of oxides and adsorbed water. The physical 
and chemical nature of these films influence the interactions between the lubricant and 
the solid surface. 

Lubricating fluids, such as mineral oils, contain small amounts of materials, such 
as fatty acids, which have attached polar (head) and hydrocarbon (tail) groups. Some 
hydrocarbon groups may contain chemically active atoms like sulfur, oxygen, and nitro- 
gen. Also, mineral oils will contain small amovints of dissolved air and water. Addi- 
tive oils may contain a variety of other elements which can participate in chemical reac- 
tions. These active components of lubricating fluids and their interactions with their 
environment determine lubrication performance. 

An ideal mechanism for lubrication by adsorbed films is shown in figure 80, where 
the interface between two opposing films provides a low-friction shear plane between the 
sliding surfaces. The actual lubricating film, however, may not be made up entirely of 
the surface active material but may be a mixture of a surfactant (e. g. , fatty acid) and a 
hydrocarbon of the base fluid. Evidence of the influence the base fluid solvent has on the 
performance of adsorbed boundary films is shown in figure 81. The lubricating ability 
is greatest when the polar component and the hydrocarbon carrier have equal chain 
lengths. 

High temperatures from the heat generated by plastic deformation and fluid shear 
cause the boundary film to desorb, disorient, or melt. This can lead to local adhesion 
between solid surfaces, as shown in figure 82. Lubricating ability is therefore fre- 
quently measured by the temperature at which the coefficient of friction rises (transition 
temperature). In general, physically and chemically adsorbed boundary films are effec- 
tive primarily at low and moderate temperatures. 

Stick-slip friction . - Stick slip is a phenomenon of friction which is characterized 
by "jerky" motion and "squealing" noise as one surface is moving relative to the other. 
Although this phenomenon is not fully understood, it is believed to be a result of small, 
local, dry, asperity contacts which are subsequently sheared. 

Fundamental Wear Concepts 

Wear, like friction, can be very complex; however, for the purpose of this discus- 
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Sion, it can be described in a relatively straightforward manner. Perhaps the most use- 
ful survey of wear mechanisms was published by Burwell (ref. 51) although other per- 
tinent publications are typified by that of Archard and Hirst (ref. 52). The various types 
of wear are frequently included in the self-description categories (table IV) of adhesive, 
corrosive, abrasive, fatigue, erosion, impact chipping, and miscellaneous (e. g. , elec- 
trical attrition). Usually, adhesive, corrosive, and abrasive wear phenomena are of 
predominant importance. However, in some applications a seal can experience all these 
wear modes. Adhesive wear phenomena are particularly dependent on the same factors 
important to the adhesive concept of friction. This should not necessarily infer, how- 
ever, that we should expect a direct correlation between wear and friction. Very fre- 
quently there is no correlation whatever. The adhesive and abrasive wear mechanisms 
are illustrated in figure 83. 

Difficulties in maintaining an effective lubricating film suggest that adhesive wear is 
of primary importance. The high ambient and higher interface temperatures achieved 
in the presence of active materials from either the lubricant or the operating environ- 
ment also suggest that corrosive wear can present a major problem. Since wear mech- 
anisms are not well understood, wear models are not usuaUy sufficiently general to 
predict quantitative values. 

Wear can be correlated on the basis of the simple theory of adhesive wear as sug- 
gested by Burwell, Archard, and others. In this model, 

V = c^WL 


and 


3H 

where V is the volume of material removed, W is the applied load, L is the total 
rubbing distance, c is the specific wear rate, C is the wear coefficient, and H is 
the indentation hardness of the softer material in contact. The wear coefficients are 
found empirically. 

Thus, the wear rates are proportional to load and inversely proportional to the 
hardness or yield stress of the worn material. Generally, wear is more severe in the 
softer material. Also wear usually increases with temperature. Some speculate that 
this is due to the decrease in hardness with temperature. Surface contamination also 
influences the wear coefficients (fig. 84). 

Rubbing contact - ultra-thin-film operation . - A broad group of seals operate under 
non-full-fibn conditions and fall into the classifications of boundary-lubricated, thin 
(mixed)-film-lubricated, or dry seals. Operation under these conditions does not lend 
itself to the straightforward treatment possible with full-film seals. This is because 
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boundary and thin-film lubrication phenomena involve the physical-chemical properties 
of the lubricant and the sliding parts. These properties are far more difficult to treat 
than the purely physical properties considered in full-film lubrication. 

Pressure-velocity criteria . - In the absence of any simple accurate analysis, 
operation under boundary lubricating conditions is frequently analyzed by the P^V (pres- 
sure times velocity) or rubbing -factor method. This method is also used as a perform- 
ance limit measure for fully lubricated face-seal operation. Although this treatment 
leaves much to be desired, it is useful if its limitations are kept in mind. 

This rubbing -fact or analysis generally is not intended to answer the question of how 
a seal will perform at any given operating condition. It is intended primarily to provide 
the limitii^ conditions under which the seal will operate satisfactorily. 

The rubbing factor is useful because it is related to two important criteria which 
can be used to establish limits for acceptable seal performance. These are the allow- 
able operating temperature and the allowable wear rate. 

Pressure-velocity-limiting criterion based on heat generation: If a reference tem- 
perature such as the ambient temperature is taken as a base, the allowable operating 
temperature can be expressed in terms of the allowable temperature rise above am- 
bient. This highly simplified approach for temperature rise calculation is aided by 
referring to figure 85. 

The heat generated under steady conditions is given by 


Qc 


F V 
= _? 2 _ 

J 


Assume heat dissipation depends on the interface temperature rise above ambient; 


= kA(T - T^) = kA AT 


Balance the heat generated with the heat dissipated and solve for the temperature rise: 

AT = — 3- 
JkA 


Note 


F 


_J2 

A 
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Thus, 


??PfV 
AT = — — 

Jk 

Letting AT be the maximum allowable temperature rise and assuming a constant coef- 
ficient of friction tj yield 


P V = ilLAI = Constant 
V 


Pressure-velocity-limiting criterion based on wear: From the relation for adhesive 
wear, the voliune of material worn from the contactii^ simface in figure 86 is given by 




Since is related to the hardness by the relation 


Real area of contact = A = — — 

° H 


H 


Dividing by A and 


and notii^ that V = Ah' and L/t = V yield 


V _ Ah' 
At At 


CwWL C^PfV 


HAt 


H 


K h'/t is the maximiun allowable wear rate specified for a particular design and 
and H are assumed to be constant, 




h'H 

tc 

w 


= Constant 


Generally, the constant based on allowable wear rate will have a lower value than the 
correspondii^ limiting constant based on heat generation. This is generally true for 
carbon seals. 
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Pressure-velocity values for various materials: Although manufacturers of bearing 
and seal materials rate their products in terms of pressure times velocity PV (usually 
in psi and ft/min), the various PV values are difficult to compare because the exact 
conditions under which they were obtained are usually not given. The advertised PV 
ratings of a wide range of materials are plotted against the heat-generation parameter 
k AT/tj in figure 87. This figure is taken from reference 41. 

Permissible values of P^V may vary with boundary-lubricated, mixed-, and full- 
film operation. It is important to check and be assured that these values have not been 
exceeded in an application. 

Application of pressure-velocity criteria: Applyir^ the PV approach to a seal 
usually means that the velocity can be replaced by its equivalent in terms of rotational 
speed V = ffDN and the temperature rise can be expressed as 


AT 


itDtjP^N 

Jk 


If k is assumed to be a constant, the temperature rise of a given-diameter seal will 
be directly proportional to rj and P^N. This relation is graphically illustrated in fig- 
ure 88. The larger the friction coefficient, the steeper the surface and the greater the 
temperatime rise for a given value of P^N (or PV). The maximiun temperature rise 
depends on the material. If the limit is exceeded, surface and internal structural dam- 
age may result. This is especially true for elastomeric lip seals. 

If the allowable operating temperature is assiuned to be fixed, the allowable temper- 
ature rise above ambient will be constant and is represented by the line AB. The pro- 
jection of AB on the base gives CD, which is the PjN = Constant in Pj and N coor- 
dinates. 


SEAL MATERIALS 

There is a great deal of information available on seal materials from both the 
technical literature and seal material suppliers. The actual material choice depends on 
the specific seal type and environment. (Generally, seal t 3 q)es are characterized as 
contact and clearance seals. ) In selecting seal materials such factors as the fluid 
cleanliness (infiltration of particles, abrasives), the temperature and chemistry of the 
fluid environment, the lubricating properties of the materials, the lubricating ability of 
the sealed fluid, and the expected or necessary life of the seal must be considered. An 
excellent survey of seal materials for specific contact- and clearance-type seals oper- 
ating in specific environments can be found in reference 53. 
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Carbon and carbon-graphite combinations are widely used seal materials. Carbon 
is found in two forms - amorphous and graphitic. Graphitic carbon has excellent bound- 
ary lubricating properties (self lubricating) but relatively low strength (low friction but 
high wear at h^h loads, and poor abrasion resistance). Amorphous carbon has high 
strength and hardness but is only a fair lubricant (low wear but higher friction). Hence, 
a combination of amorphous and graphitic carbon is a useful compromise of self- 
lub 'icating properties and relatively high compressive strength. The coefficient of 
friction at room temperature is about 0. 10 to 0. 75, and comprehensive strengths are 
103 to 276 MN (15 000 to 40 000 psi). Additional useful properties are low density, low 
elastic modulus (conformability), and resistance to thermal shock and corrosion. Most 
physical properties of carbon improve with temperature. Ejq)erimenters have observed 
that oxygen and/ or water vapor play an important role in the physical behavior of carbon 
and graphite. The absorbed gases act to reduce friction between the crystallites of the 
carbon body and those crystallites of the carbon adhering to the matii^ surface. It is 
this ability to develop transfer films on opposing materials that gives carbon its excel- 
lent friction and wear properties. The overall conclusion (generalizing results) is that 
the matii^ material and its surface finish play a major role in the frictional behavior of 
a carbon. Also load level, speed, and ability to dissipate heat from the mating surface 
are important. 

Since the carbon surface is porous, it can be impregnated with oxidation inhibitors 
(oxidation being a problem at high temperatures) or treated to prevent air from permeat- 
ing into the carbon. Treatments are also used to improve the friction and wear proper- 
ties and the chemical resistivity. 

Because carbon has a low modulus of elasticity (about 27 580 MN (4x10® psi)), it 
readily conforms to the mating surface. This property, aloi^ with its excellent friction 
and wear properties, has led some people to call carbon a very "forgiving" material. 
More information on carbon and graphite is given in reference 54. 


CONCLUDING REMARKS 

Some fvindamentals of fluid sealing have been presented. General fluid-flow equa- 
tions were developed that are applicable for many fluid seals. Included in this report 
are seal performance parameters and seal operating regimes. Also discussed are ge- 
ometry effects, surface deformations, rotation, pressure balancing, mechanisms of 
seal surface separation, fundamental friction and wear concepts applicable to seals, 
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seal materials, and pressure-velocity criteria. These fundamentals should enhance 
imderstandir^ of various sealii^ devices and enable more in-depth investigation of spe- 
cific seals. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 2, 1975, 

505-04. 
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APPENDIX A 


SYMBOLS 

A cross-sectional area or apparent area of contact, cm^; in. ^ 

A^ real area of contact 

a speed of sound, m/sec; ft/sec 

B friction parameter, 4fL/D 

C conversion constant 

Cj) flow discharge coefficient 

Cj^ velocity entrance loss coefficient 

wear coefficient 

Cp specific heat at constant pressure, J/(kg)(K); Btu/(lbm)(°R) 

c^ specific heat at constant volume, J/(1^)(K); Btu/(lbm)(°R) 

c^ specific wear rate, cm^cm-N; in. Vin- -Ibf 
D hydraulic diameter, 2h 

D/Dt material derivative, d/dt + u(a/3x) + v(a/3z) + w(3/ay) 

d diameter, cm; in. 

e eccentricity; or amplitude of normal movement, cm; in. 

e' specific internal energy 

F seal opening force, N; Ibf 

F pressure profile load factor, dimensionless 

f Fanning friction factor, T^/(pu^/2) 

f mean Fanning friction factor 

G flow parameter, M/A, kg/min/m^; Ibm/sec/m^ 

g acceleration of gravity, m/sec^; ft/sec^ 

H indentation hardness; or total head loss coefficient 

HP horsepower 

h film thickness (gap), cm; in. 

Ah change in film thickness, pm, mil 

h’ mean depth of surface worn away, cm; in. 
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^char 

^n 

i 

J 

K 

Kn 

k 

k’ 

L 

L 

^2 

M 

M 

AM 

m 

N 

N 

P 

AP 

P 

^min 

Pr 

Q 

% 

% 

q 

R 


characteristic film thickness, ^ 

minimum film thickness, cm; in. 
specific enthalpy 

heat equivalent of mechanical work conversion factor, 4. 184 J/cal (778 ft-lb/ 
Btu) 

expression defined by eq. (157a); or axial film stiffness 
Knudsen number, xA 

thermal conductivity; or Boltzmann constant 
head loss coefficient 

flow length from entrance to point of choking 
total rubbing distance, cm; in. 
length of recess region in step bearing 
length of land region in step bearing 
Mach number 

mass leakage flow rate, kg/min; Ibm/min 
change in mass flow 

molecular weight of gas, kg/kg-mole; Ibm/lbm-mole 

rotational speed, rpm; or number of flow restrictors 

molecular flow rate, molecules/sec 

Avogadro number, 6. 027x10^^ molecules/mole 

static pressure, N/m ; psi 

pressure differential, N/m^; psi 

yield strength of material, N/m^; psi 

smaller pressure of two pressure boundary conditions, N/m^; psi 
normal force, N; Ibf 

volximetric leakage flow rate, scms; scfm 
heat dissipated, J/sec; Btu/sec 
heat generated, J/sec; Btu/sec 

volumetric flow rate per unit width, m^/sec/m; ft^/min/ft 
radius, cm; in. 
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91 gas constant, ^m, m-N/(kg)(K); in. -lbf/(lbm)(°R) 

R universal gas constant, 8. 3143 N-m/(mole)(K); 1545. 4 ft-lbf/(lbm-mole)(°R) 

R mean radius, (Rj + R2)/2, cm; in. 

AR seal face length 

Re leakage-flow Reynolds number in radial direction, pUh/jn 

Re^. shear-flow Reynolds number in circumferential direction, pRS2h/ju 

r radial direction coordinate 

o 

S shear strength of material, N/m ; psi 
s labyrinth knife-edge spacing, cm; in. 

T temperature, K; °F 

^ linear tilt factor, h 2 ( 2 hj + 

AT temperature rise in fluid film, K; °F 
Ta Taylor number 

t labyrinth knife-edge thickness, cm; in. ; or time, sec 

U leakage-flow reference velocity, m/ sec; ft/sec 

Uj lower moving surface velocity, m/sec; ft/sec 

U 2 upper moving surface velocity, m/sec; ft/sec 

u velocity in r-direction or x-direction, m/sec; ft/sec 

V reference flow velocity across gap, U(h/AR); or seal rubbing velocity, m/sec; 

ft/ sec 

— 3 3 

V volume of wear, cm ; in. 

V velocity in z-direction, m/sec; ft/sec 

W flow width, cm; in. 

W normal load, N; Ibf 

W reference shear -flow velocity (surface speed due to rotation), m/sec; ft/sec 

w velocity in 0 -direction or y-direction, m/sec; ft/sec 

3 3 

X body force in x-direction, N/cm ; Ibf/ in. 

Xg center of pressure in radial or x-direction, cm; in. 

X coordinate in pressure gradient direction 

Xg entrance length, cm; in. 
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V body force in y-direction, N/cm^; Ibf/cm^; or expansion function 
y distance normal to flow, cm; in. 

O O 

Z body force in z-direction, N/cm^; Ibf/in. 
a seal-face distortion (relative inclination angle of surface), rad 
a fluid acceleration 

/3 radial flow factor, W/r; or ratio of nozzle-inlet- to nozzle -exit-flow cross- 

sectional area in section INVISCID FLOW EQUATIONS; or ratio of exit to inlet 
film thickness in section Flow Between Flat Convergii^ or Diverging Plates, 

Vi 

y specific-heat ratio, Cp/c^ 

6 bovmdary layer thickness, cm; in. 

e excursion ratio, e/hj^^ 

rj coefficient of Coulomb friction 

0 circumferential coordinate 

A bearing or compressibility number, 6/xU 

X molecular mean free path, cm; in. 

fi absolute or dynamic viscosity, N-sec/cm ; Ibf-sec/in. 

V kinematic viscosity, m /sec; ft /sec 

p density, N-sec^cm^; Ibf-sec/in. ^ 

CT squeeze number, 12 

T shear stress, N/cm^; psi; or characteristic dimensionless time 

(p ratio of inlet to exit film thickness 

0 angular rotational velocity, rad/ sec 

w normal relative surface oscillation frequency, Hz 
2 

V Laplacian operator 
Subscripts: 

a ambient conditions 

av average 

c cavitation or degasification 

char characteristic 

F frictional 
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f net face contact 

HS hydrostatic 

h based on film thickness 

I inertial 

inc incompressible 

L entrance length 

m mean 

max maximiun 

min minimum 

r based on rotational velocity 

ref reference condition 

S spring 

SD sealing dam 

s step location 

std standard conditions 

w wetted surface 

X x-direction 

z z -direction 

7 ] shear force (friction) component 

0 sealed (reservoir) conditions 

1 entrance conditions 

2 exit conditions 

3 ambient sump conditions 
Superscript: 

* dimensionless quantity; or reference to Mach 1 condition (critical flow) 



APPENDIX B 


VISCOSITY OF FLUIDS 

Viscosity is the property of a substance by which it offers resistance to shear 
stresses. In a Newtonian fluid, the »riscosity is in linear proportion to the ability of the 
fluid to resist such shear stresses. 

All substances, both liquids and gases, have viscosity. Because of this, the prop- 
erty should be explained by one or more physical properties common to all fluids, such 
as molecular activity. Because of molecular activity, there is a constant interchar^e 
of molecules, and therefore of momentum, between adjacent layers of the fluid. If ad- 
jacent layers are moving with different velocities, this constant interchange of momen- 
tum sets up a resistance to any relative motion of the two layers (fig. 89). As a result 
of this resistance, energy is transformed into heat; and to maintain the velocity, a 
steady force is required. 

The viscosity of a liquid decreases with an increase in temperature. Unlike other 
material properties which change slightly with temperature, viscosity is somewhat 
unique in that, for a liquid, it can change several orders of magnitude with temperature. 
The temperature range seen by a liquid in a hydrodynamic environment is usually very 
large; thus, predicting the variations of viscosity with temperature is important. The 
interchai^e of momentum is accelerated with an increase in temperature, but the vis- 
cosity of a liquid must be regarded as the combined effect of cohesion and interchange of 
momentum. 

Cohesion is the force with which like molecules of a substance attract each other. 

A change in temperature has opposite effects on cohesion and molecular activity, with 
the effect of cohesion being more pronounced so that in a liquid as temperature is de- 
creased, viscosity is increased. Because of low cohesive forces in gases, the viscosity 
of a gas increases with temperature. The viscosities of some liquids and gases are 
compared in figure 90. 

Note the analogy with solids, where Hooke's law relates normal stress with normal 
strain; for example, Young's modulus is a material property. For fluids the stress is 
proportional to the strain rate. (When the fluid is not in motion, the shear stresses 
vanish. ) 

The simple linear dependence of shear stress on velocity gradient is known as the 
Newtonian relation. Fluids which behave this way are called Newtonian fluids. Many 
oils at moderate shearing rates, water, air, and many gases are essentially Newtonian. 
Thus, the theory is highly developed for these fluids. (All equations and formulas in this 
text are for Newtonian fluids. ) 

Seals may use fluids that are not Newtonian. Some fluids are a function not only of 
velocity gradient but of ordinary strain as well. These fluids are termed "viscoelas- 
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tic. " (Polymer-fortified oils are viscoelastic, especially at high shear rates. ) Some 
fluids are not Newtonian because there may be a rather complicated nonlinear shear 
stress - strain rate relation. Some greases and sludges behave as solids until an ap- 
parent yield stress is exceeded; then they behave as viscous fluids. This is known as 
plastic behavior. Various shear stress - strain rate relations are depicted in figure 91. 
Further information on non-Newtonian behavior in seals can be found in reference 2. 
Understanding of Newtonian behavior can shed light on the physical behavior of non- 
Newtonian fluids. 
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APPENDIX C 


COMPRESSIBLE FLOW - AN ALTERNATE THEORETICAL APPROACH TO 
FINDING THE ORIFICE EXPANSION FUNCTION 


Consider the adiabatic one -dimensional energy conservation equation. (Again for 
compressible flow the energy equation must be considered. ) 

P u2 

e' + — + — = Constant 
P 2 

or 

2 2 
u. u„ 

ii + — = i.) + — = Constant (Cl) 

^2 ^2 


For an ideal gas, 

i = CpT 

And by neglecting Uj, equation (Cl) becomes 

U2 = ^2Cp(Ti - Tg) 

Now 


c 

P y.l 


Thus, 


u 


1 ■ 




Since 



the exit velocity is 





Y = 


Thus, the orifice expansion function is found theoretically. However, in practice, it has 
been found to generally underpredict the flow rate when a constant discharge coefficient 
is assumed. 
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APPENDIX D 


DERIVATION OF REYNOLDS LUBRICATION EQUATION^ 


As stated in the main text, load capacity can be built up in a hydrodynamic bearing 
from both the we(%e-film and the squeeze-film actions. 

The differential equation that describes hydrodynamic bearing operation relates 
viscosity, film thickness, pressure, and sliding velocity. This equation was first de- 
veloped by Osborne Reynolds and is named for him. In the derivation of the Reynolds 
equation, the following assumptions (also used in deriving the equations in the section 
Mechanism of Film Pressure Generation Between Relatively Moving Surfaces) were 
used: 

(1) The flow is laminar. 

(2) The fluid is Newtonian (Newton's definition of viscosity. Viscosity = Shear 
stress/Shear rate, holds), and the viscosity is constant across the film, 

(3) Inertial forces are small as compared with viscous shear forces and are ne- 
glected, (Reynolds number is low. ) 

(4) The weight of the fluid is negligible compared with other forces. 

(5) No slip exists at the fluid- solid interface. 

(6) The film thickness is small compared with any radii of curvature. 

(7) The pressure across the film is constant. 

(8) The inclination of one surface relative to the other is so small that the sine of the 
angle of inclination can be set equal to the angle and the cosine can be set equal to unity. 
Figure 92 shows that the moving member has a tangential velocity U and a normal ve- 
locity V. 

The coordinate axes are directed as shown, with the origin in the moving member. 
An infinitesimal volume of fluid, with dimensions dx, dy, and dz in the film, is shown. 
Pressure forces act on four faces and shear forces act on two faces. There is no net 
pressure force in the y-direction (assumption 7). Let the fluid particle have velocities 
u, V, and w in the x-, y-, and z -directions. 

If the forces in the x-direction are equated with negligible fluid accelerations. 


9t, 


3P 


P dy dz + (T + — - dy ) dx dz - ( P + — dx) dy dz - t dx dz = 0 
\^dy \dx 


(Dl) 



9x 9y 


^Appendix is partially derived from ref. 40. 


(D2) 
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From the definition of viscosity, 


Similarly, in the z-direction 


3u 

9y 


^ 3y2 


ap , a^w 

— = M — - 

0y2 


(D3) 


(D4) 


(D5) 


E^quations (D4) and (D5) can be integrated with respect to y since pressure P is inde- 
pendent of y. Ihtegratii^ twice yields 


u = Cl + C2y + — y^ 

^ 2p ax 

(D6) 

w - Cg + C4y + y2 

^ 2p az 

(D7) 

Specific boxmdary conditions are used to determine the constants 

Cj, Cg, Cg, and 

C4. Li the x-direction, Ci and C2 are determined by 

u = U at y = 0 


u = 0 at y = h 


where h is the thickness of the film at any point. Then 


u(y) - - y) u - ~ y) 

(D8) 

h 2/i 3x 


The variation of the velocity u across the film as given in equation (D8) is the sum 
of two terms. The first is the linear variation in velocity from the stationary surface at 
y = h to the surface moving with velocity U at y = 0. It is independent of viscosity and 
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is constant throughout the film. The second term in the velocity distribution is that due 
to the pressure gradient in the x-direction. The distribution across the film is parabolic 
and is dependent on viscosity. 

In the z-direction, Cg and are determined by 

w = 0 at y = 0 

w = 0 at y = h 


Then 


w(y)= (D9) 

2p 3z 

In the z-direction, the velocity distribution is parabolic and results from the pressure 
gradient in that direction. 

Because the Reynolds equation is essentially a continuity equation, the next step in 
its derivation is the calculation of flow rates in the film. When a unit width is consid- 
ered, the flows in the x- and the z-directions, q^ and q^, are obtained by integration of 
the velocities u and w over the film thickness: 



h^ 8P 
12 n dx 


(DIO) 


^z 



dy = 


h^ 3P 
12ju 0z 


(Dll) 


If a volume of fluid of cross section dx dz extends across the film of thickness h and 
an expression is written for the net mass flow out of the volume in the x- and z- 
directions, this expression must equal the rate of chaise of mass in the volume (fig. 93): 

!5^dxdz.^dxdz = -iMdxdz (D12) 

9x dz 3t 


where 8h/9t = V is positive when away from the other surface. Substituting for q^ 
and q„ from equations (DIO) and (Dll) gives 
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Equation (D13) is a general form of the Reynolds equations. It applies for full-film 
lubrication in many bearing configurations for both compressible and incompressible 
fluids when the assumptions made in its derivation are valid. Solution of the Reynolds 
equation for a particular bearing configuration results in an expression for the pressure 
throughout the film as a function of viscosity, film shape, and velocity. Integration of 
the pressure over the film area yields the total load supported by the film. 

No exact solution of the generalized form of the Reynolds equation given by equa- 
tion (D13) has been obtained. Approximate solutions for specific conditions have been 
obtained by the use of infinite series and relaxation techniques, and accurate numerical 
solutions have been obtained through the use of digital computers. While all these ap- 
proximate solutions are valuable, they fall far short of a general solution. 

The Reynolds equation is readily solvable if it is assumed that the bearing is wide 
relative to its length (the so-called infinite bearing) or that the bearing is short enough 
that the pressure-induced circumferential flow is negligible (the short-bearing solution). 
Dowson (ref. 56) has developed a generalized Reynolds equation for fluid-film lubrica- 
tion. 
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TABLE I. - OTHER NAMES FOUND IN THE 
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‘Ref. 6 also defines a complement to unity which 
is called the degree of balance. 






TABLE m. - MASS LEAKAGE FLOW RATE EQUATIONS FOR VARIOUS QUASI- FULLY- 


DEVELOPED (SUBSONIC) FLOW RATE SITUATIONS 


Case 

Laminar flow 

Turbulent flow 

Constant area, parallel 
surfaces 

wh3(p;-Pi) 

Wh^^'’^(^Pj - p|)"*' 

JVl - ' — 

24/i£»T(R2 - R^) 

2^ ' '^(0. 079)‘^ ''^«‘*'''^T^^'^;i^'^'^(R2 - Rp'^' 

Radial flow 


34/7,12/7(p2 _ p2)^/'^ 

12p^T In — 
^2 

2(0.079)4/^1/7^4/7^4/7/ 1 . 1 ] 

h3/4; 

Small linear deformations 
(h = + ax) and constant 

width 


W.;V7(p2 . ’ 

24M£fT(R2 - R^) 

2®/1(0. 079)4/1^ l/1flj4/7.p4/7^p^ _ R^4/7 


TABLE IV. - TYPES OF WEAR (FROM REF. 40) 


Type of wear 

Example 

Adhesive 

Interface welds 

Corrosive 

Chemical surface reactions 

Abrasive 

Rough surface and free-particle cutting 

Fatigue 

Subsurface repetitive-shear-stress failure 

Erosion 

Impact displacement by fluid (possible with 


entrained solids) 

Impact chipping 

Solid-surface-impact stress failure 

Miscellaneous 

Electrical attrition and other phenomena 
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Low load 



Film thickness 
same order 
as surface 
roughness 


High load 



Laminar flow 
disturbed 


Very high load 



Figure 2. - Hydrodynamic, thin-film, and boundary lubrication modes. 
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Compressibility 

effects 


Turbulence- 


Seal duty parameter, 

Figure 3. - Variation of friction coefficient for gas lubrication, compared with liquid 
lubrication. 




• Laminar flow * 




Figure 4. - Stressed element in fluid film. (From ref. 5. ) 





Secondary 
seal diameter 



Figure 7. - Pressure -balanced face seal (no axial film stiffness). 



Figure 8. - Radial face seal, identifying nomenclature. 




2)U dx 



A B C (simple D 


Couette flow) 



> 0 

dP 

> 0 

o 

II 

^<0 

dx 


dx 


dx 

dx 

du 

< 0 

du 

= 0 

^ >0 

du I 

dz 

!z=o 

dz 

z=0 

‘^z=0 

o 

11 

N 


Figure 10. - General Couette flow velocity profiles, where 

u = J-^(z2-hz)+U5. 

2/1 dx h 
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Figure 13. - Representative pressure profiles for var- 
ious surface conditions, where p is ratio of exit to 
inlet film thickness, 


P 



Figure 14. - Linear converging surfaces of a simple slider 
bearing in relative motion, illustrating pressure distri- 
bution. 


Ill 





T T + dT 
p P + dP 
M M + dM 


Figure 17. - Model of seal faces, showing notation and including control volume for 
quasi“One-dimensional flow with area change. 
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Figure 18. “ Effect of eccentricity on turbulent flow. (From 
ref. 10). 



Figure 19. - Comparison of effect 
of eccentricity on laminar and 
turbulent flow. (From ref. 10, ) 






(a-2) Diverging (Ah = - hj < 0; = hg). 

(a) Linear face contours. 



(b) Parabolic face contours. 


Figure 20. - Seal -face contours studied by Snapp (ref. 14). (Analytical 
film thickness expressions are given in ref. 14. ) 
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0 600 1200 1800 2400 3000 

Inlet water pressure, Pp psi 


Figure 21. - Minimum film thickness as function 
of inlet water pressure for various face con- 
tours. Leakage flow rate, Q, 0.0019 m /min 
{0. 5 gal/min); change in film thickness, Ah, 
12.7 jam (0. 0005 in. ). (From ref. 14; re- 
vised. ) 
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Face contour 
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Ratio of leakage flow rate to inlet water pressure, 
Q/Pj, (m^/min)/(N/m^) 


Ratio of leakage flow rate to inlet water pressure, 
Q/Pp (gal/min)/1000 psi) 


Figure 22. - Minimum film thickness as function of ratio of leakage flow rate to inlet 
water pressure for various face contours. (From ref. 14; revised.) 
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Film pressure, psi 

Figure 23. - Radial pressure profiles for various face con«- 
tours. (From ref. 14; revised. ) 


118 



\\\ 


V > rt 


N[ ‘M "P^oi TBUiJONT 


00 o CD eg 00 

CM eg C^J 1-H rH 


w \ 

u bii \ 

^ .S w) Xv V 

s & -2 xO 


X\ 




jqi M pBoi iBuuoN 


119 


vised. ) 



Face contour 





o 


^ t > 

^ c 
a o 


1 I 3 

3 “ 

:=5 ^ a) c« 3 

rt -^3 > O a; 

^ O S S G 

rt O o .3 

A O O J 


. XL 
o 




c 

S 



C? 



0> 

ij 


M "M ‘P^OT TBUIJON 




J ^ ^ I 


jqi ‘M ‘PTOX lEOTJON 


120 


Figure 26. - Seal face load as function of leakage flow rate 
for various face contours. Constant exit pressure; con- 
stant film thickness. (From ref. 14; revised.) 
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Figure 29. - Load factor as function 
of pressure ratio in subsonic flow 
regime for parallel sealing sur- 
faces. 



Figure 30. - Density and velocity as 
function of flow distance. 



Air leakage flow rate, M, Ibm/sec 



Air pressure, psig 


Figure 31. - Labyrinth seal leakage as function of air pressure. (From 
ref. 16.) 


Seal-face radial length, AR, mil 


Sealed-gas 

pressure, 

■ Po> 

p 
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Compressible flow 
Incompressible flow 



Sealed-gas 
temperature, 
K(°F) / 
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-// 


148 (215) , 


45 (65) 


Pressure ratio, P/P 3 

32, - Pressure profiles of gas within parallel seal faces. Sump pressure, P 
2 

M/cm abs (15 psia); subsonic flow (Mach number < 0, 845), 








Seal-face radial length, mil 



Pressure ratio, P/P 3 


Figure 35. • Pressure ratio as function of mean film thickness. Sealed-gas pressure, Pq, 

45 N/cm2 abs (65 psla); sealed-gas temperature, Tq, 311 K (100° F); sump pressure, P 3 , 
10.3 N/cm^ abs (15 psia); seal-face radial length, AR, 0. 127 cm (50 mils). 
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Figure 37. - Comparison of classical viscous flow theory and quasi-one-dimensional analysis with ra- 
dial flow experimental results. Sealed-gas pressure, Pq, 41. 8 N/cm^ abs (60. 6 psia); sealed-gas 
temperature, Tq, 300 K (80° F); sump pressure, Pg, varied; film thickness, h, 38 ^m (1. 5 mils); 
entrance radius, Rj, 6.98 cm (2.75 in.); exit radius, R 2 , 7.62 cm (3 in.). 
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Mass leakage flow rate, M, Ibm/min 



Figure 39. - Mass leakage flow rate as function of pressure ratio for vari- 
ous seal-face radial lengths. Sealed-gas pressure, Pq, 148 N/ cm abs 
(215 psia); sealed-gas temperature, Tq, 700 K (800*^ F); film thickness, 
h, 10. 2 pm (0. 4 mil); entrance radius, Rj, 8. 293 cm (3. 265 in. ). 
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Figure 42. - Seal opening force as function of film thickness for parallel faces; comparison of classi- 
cal viscous flow theory and quasi-one-dimensional analysis. Sealed-gas pressure, Pq, 

148 N/cm^ abs (215 psia); sealed-gas temperature, Tq, 700 K (800^ F); sump pressure, P^, 

10.3 N/cm^ abs (15 psia); seal-face radial length, AR, 0, 127 cm (50 mils). 
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Discharge coefficient, 



Pressure ratio, P/ Pq 


Figure 45. - Discharge coefficient as function of pressure ratio for nozzles and ori 
fices. (From ref. 20; revised. ) 



Rounded 

orifice 



(a) Low- velocity flow through round- 
edge orifice. 



(b) Low- velocity flow through sharp-edge 
orifice. 



(c) High-velocity flow through sharp- edge 
orifice. 


Figure 46. - Aerodynamic flow characteristics of air. (From ref. 21.) 



Figure 47. - Schematic representation of orifice 
flow driven by gravitational head. 



Figure 48, ~ Schematic of typical subsonic nozzle. 








(a) Laminar flow in entry region of pipe. 



Figure 50. - Qualitative comparison of entrance lengths of laminar and tur- 
bulent flows with smooth (rounded) entrance lengths. 



Figure 51. - Entrance loss mechanisms: turning (separated flow); wall 
shear (boundary layer growth); and momentum increase (less uniform 
velocity distribution). 
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Dimensionless velocity, U/U 


Center 



Dimensionless distance, z/(h/2) 


Figure 54. - Qualitative comparison of laminar 
and turbulent velocity distribution. 



Flow velocity, V 


Figure 55. -Representation of pipe pressure drop for laminar, transition, 
and turbulent flow. 




Radial velocity function Dimensionless circumfer- 

ential velocity 



(a) Circumferential ve- 


locity profiles. 



across gap, z* 

(b) Radial velocity pro- 
files. 

Figure 58. - Theoretical 
radial and circumfer- 
ential velocity profiles 
for various Reynolds 
numbers. (From 
ref. 33.) 
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Figure 59. - Comparison of theoreti- 
cal and experimental velocity pro- 
files for boundary layer case 
(Re = 441). (From ref. 33.) 





streamline 

direction 


K \ V 


(b) Closeup of vortex flow, with diagram indicating 
streamline directions. 

Figure 61. - Flow visualization of fully developed Taylor 
vortices for concentric cylinders with no axial flow. 
(From ref. 34. ) 




Rotational Reynolds number (circumfer- Mean axial velocity, 



Figure 62. - Schematic representation of different flow re- 
gimes in an annulus with axial flow and inner cylinder ro- 
tating. (From ref. 36.) 
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X//////A Quasi-incompressible flow 
1 I Transition (laminar to turbulent 

flow) 

Turbulent flow (radial and shear 
flow coupled) 


^900 



' W Theoretical limit 
of isothermal flow 

t 7“ Choked flow 
/ (laminar or 
t t urbu lent ) 


1/Vy 1 

Mach number (hydrostatic radial 
pressure flow direction), M 


Figure 63. - Envelope of possible seal-face flow regimes 
for gases. (From ref. 38.) 
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(a) Hydrodynamic. 


Recess ^ Fluid ex- 

n pelled 

f ^-Pressurized flow 
from external source 
(b) Hydrostatic. 
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I fluid pressure 
distribution 
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Fluid ex- 
pelled 
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pelled 



(c) Squeeze film. 


- Pressurized flow 
from external source 
(d) Hybrid, 


Figure 64. - Modes of fluid-film lubrication. (From ref. 39). 
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Figure 65, - Velocity distribution between two parallel 
plane surfaces in relative motion. (From ref. 40. ) 
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u 



(a) Velocity distribution as a result of relative motion. 



(b) Velocity distribution as a result of pressure 
buildup in lubricant film. 

^.0 



(c) Velocity distribution as a result of superimposed 
velocity and pressure flows. 

Figure 66. - Flow between two inclined plane surfaces. 
(From ref. 40. ). 




P - 0 


(b) Step bearing pressure distribution. 
Figure 67. - Step bearing geometry and pressure distribution. 
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inc inc 




(b) Effect on friction, center of (c) Effect on load capacity, 

pressure, and load capacity, 

Xinc = 0.582 AR; 

A/AR; W.„„ = 0.069 
/iU(AR)2L/h|. 

Figure 68, - Typical compressibility effects in hydrodynamic bearings. Slider bear- 
ing; ratio of flow length to seal-face radial length, L/ AR, 1; ratio of entrance to 
exit film thickness, h^/h2, 2. (From ref. 41.) 
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Figure 69. - Flow between parallel plane surfaces that approach each 
other. (From ref. 40. ) 



0 T 2r 


cot 

Figure 70, - Dimensionless squeeze-film 
force (load as function of time) for one 
cycle, for small and large excursion 
ratios. Squeeze number, a, 10. 

(From ref. 43. ) 
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Figure 71. - Typical pressure profile found from incompres- 
sible flow solution, illustrating that "degassing” or cavita- 
tion will yield a net squeeze film force over a cycle. 



Figure 72. - Typical appearance of clearance 
between two seal surfaces, (From 
ref. 45.) 


Y Convergent sur- 
Microface -j ^ f^ce spacing, f 



Figure 73, - Idealized appearance of clearance between 
two seal surfaces, (From ref. 45. ) 
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B 


B 



Without distortion 



Figure 74. - Pressure distribution on a sym- 
metrical microface pair. (From ref. 45. ) 



Position 


Figure 75. - Pressure measurement for large-scale asperity model. Ambient 

2 2 
pressure, 177 kN/m (25. 7 psia); cavity pressure, 103 kN'/m (14. 9 psia); 

2 

peak pressure, 463 kN/m (67. 2 psia). (From ref. 46.) 
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Rotation of 
wavy face, ^ 

Figure 76. - Wavy face -seal surfaces, illustrating hydrodynamic pres 
sure generation. (From ref. 47.) 
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Figure 78. - Influence of contamination on friction coefficient for vari- 
ous material combinations. (From ref, 49. ) 



Motion 



Figure 79. - Typical solid surface in contact with lubricant - 
undisturbed surfaces. 



Motion 



Figure 80. - Lubrication by adsorbed lubricant - lower effective shear 
strei^th (lower friction than solids in contact). 



Motion 



Figure 81. - Failure of surface films - friction increases and wear occurs. 



Motion 



Extreme -pressure 
film or chemisorbed 
liquid (polymer) 


— Adsorbed lubricant 


Chemisorbed atmos- 
pheric elements 


Figure 82. - Abraded surface after contact - oxide formation rate greater 
than removal rate. 


Elastic and plastic deformation 



^ Adhesion (welding) 



Wear particle (can be abrasive) 


L Ploughing (due to 
gross plastic flow) 


Figure 83. - Adhesive and abrasive wear mechanisms. (From ref. 48.) 



Wear coefficient with 
loose particles 


Materials in contact 



high vacuum in air (water; gasoline; (pure mineral (mineral oil with 

nonwetting liquid oil; molten lubricity additive; 
metal) glass; wetting fatty oil; good 

liquid metal) synthetic lubricant) 

Figure 84, - Influence of contamination on wear coefficients for various material 
combinations. (From ref. 49. ) 
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Figure 85. - Simplified model for deriving 
the pressure-velocity-limiting criterion 
based on heat generation. 


P 



Figure 86. - Simplified model for deriving the 
pressure -velocity-limiting criterion based 
on wear. 
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(From ref. 41, ) 




Figure 88. - Graphical representation of relation 
(P^V) - Constant. (From ref. 41.) 



Figure 89. - Newton’s law of viscosity: F = ju(AU/h) or ii = (F/A)/(u/h) = Shear stress/Shear 
rate, where is viscosity, U is leakage flow reference velocity, h Is film thickness, and 
A, u, and F are defined in sketch. 
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Figure 90, - Absolute viscosity of some fluids. 



Figure 91. - Shear stress as function of shear strain rate for various 
types of viscous and plastic fluids. 
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